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INTRODUCTION 

Selection  and  ranking  (orderi nc  ■  problems  in  statistical  in  forma  .• 
arise  mainly  because  the  classical  tests  of  homogeneity  are  ntii-n  in 
adequate  in  certain  situations  where  tin-  experimenter  is  i n t « •»-» •  . t . •• !  in 
comparing  k  (•  2)  populations,  treatments  or  processes  with  the  >m.il  oi 
selecting  one  or  more  worthwhile  (oood)  populations.  MosH  t  r  (in;:;. 
Paulson  (1919),  Bahadur  11950)  and  Bahadur  and  Bobbin ^  1 1  *‘,0 )  Wnm  ,rn:i>i 

the  earliest  research  workers  to  recognize  this  inadequacy  .end  to  im 
mulate  the  problem  as  a  multiple  decision  problem  aimed  at  the  ,«■  I <m  t  im 
and  ranking  of  the  k  populations. 

In  the  thirty  years  since  these  early  papers,  selection  and  rant  inn 
problems  have  become  an  active  area  of  statistical  res.., orb .  Ihm  have 

been  two  approaches  to  these  problems,  the  'indifference  .-one'  nppMM,  p 
and  the  'subset  selection'  approach.  In  the  first  approach,  due  to 
Bechhofer  (1954),  the  experimenter  wishes  to  selec  t,  one  populal  ion 
a  fixed  number  t  •  1  of  population)  which  is  guaranteed  to  bo  the-  . im- 
of  interest  to  him  with  a  fixed  probability  P*  whenevet  t ho  urtnow. 
parameters  lie  outside  some  subspace  of  the*  parametor  spat  i*.  tin-  •... 
called  indifference  zone.  Important  contributions  u-.irui  t h i • .  appi.c.uti 
have  been  made  by  Bechhofer  and  Sobol  (1954),  Bee hhofoi  ,  buiuu-t t  ,m.i 
Sobol  (1951),  Sobol  and  lliiyctt  ( 1 <>r*  7 ) ,  Sobol  (lOfw’l,  Hoc  him  f  o>  ,  !  i.:m 


and  Sobol  (196R),  Mahaiunulu  (1%7),  D«*su  .md  Sobol  (lot,;*,  I '  >  '  1  1  ,m,| 


iamhane  (tnd  Bechhofee  (  1977  ,  1  9  79  )  among  others.  A  quite  complete 
bihl iography  may  be  found  in  Gupta  and  Panchapakesan  (1979'  (see  also 
Gibbons,  Olkin  and  Sobel  (1977)). 

The  second  approach  pioneered  by  Gupta  (1956,  1967,  1966)  assumes 
no  a  priori  information  about  the  parameter  space.  A  simile  pop  slat  i;;’' 
is  not  necessarily  chosen;  rather  a  subset  of  the  given  k  populations 
is  selected  depending  on  the  outcome  of  the  experiment.  It  is  uuaran'e 
to  (.on tain  the  population(s)  of  interest  with  probability  which  is 
least  equal  to  P*  (the  basic  probability  requirement)  regardless  of  : he 
true  unknown  configurations  of  the  parameters.  Some  recent  contributor 
in  the  category  of  subset  selection  include;  lieely  (1965),  Gnanadesiku 
(1966),  Gnanadesikan  and  Gupta  (1970),  Gupta  (1967),  Gupta  and  Studden 
(19/0),  Nagel  (1970),  Gupta  and  Nagel  (1971),  Gupta  and  Panchapakesan 
(19/P),  Rizvi  and  Sobel  (1967),  McDonald  (1969),  Gupta  and  McDonald 
(10/0),  Santner  (1975),  W.  T.  Huang  (197?),  D.  Y.  Huang  (1975),  Gupta 
and  Huang  (1975a,  1975b)  and  Gupta  and  Huang  (1976). 

Subset  selection  procedures  can  also  be  thought  of  as  screening 
procedures  which  enable  the  experimenter  to  select  a  subset  of  pop¬ 
ulations  (under  study)  which  contains  the  populations  interest  so 
that  the  populations  in  the  selected  subset  can  be  further  studies. 

Sequential  and  multistage  aspects  of  the  ranking  and  selection 
problems,  have  been  explored,  based  on  the  indifference  zmte  approach 
by  Mechhofer,  Dunne tt  and  Sobel  (1969),  liochho'er  (195K),  I’.mlson 
(196?,  |96H,  1964,  1967)  and  Bechhofer,  Kiefer  and  Sobel  (!%?’). 

Parson  and  Gupta  (197?),  Huang  (197?',  Gupta  and  Huan  ;  (19/5).  Gupt  i 
and  Miesike  (1979)  and  f.arrol  i  (19/4)  hive  investigated  s  ii  .  t  •  •  i  - 

prn< educes,  based  on  sequential  sampling. 


Contributions  to  optimum  properties  of  subset  selection  prm  educes 
have  been  made  by  Lehmann  (1961),  Studden  (1967),  Deely  and  Gupta  (Iugk) 


Berger  (1977,  1979),  Gupta  and  Hsu  (1978),  Gupta  and  MiesUo  ( I •*/;'. .) . 
Berger  and  Gupta  (1980). 

In  the  decision-theoretic  approach  to  the  subset  selection  problems, 
Goel  and  Rubin  (1977),  Chernoff  and  Yahav  (1977),  Bickel  and  Ynhuv  ( I h// 
Gupta  and  Hsu  (1978),  Miescke  (1979),  Gupta  and  Kim  (1980),  'lugia  and 
Hsiao  (1980)  have  given  different  formulations  under  different  lev. 
functions  and  carried  out  investigations  which  indicate  that  the  Gupta 
type  maximum  (minimum)  means  procedures  are  quite  'optimal1  md  ' :  id  ns  t  ' 
The  main  purpose  of  this  thesis  is  to  study  seme  problems  us  mi 
the  subset  selection  approach  and  provide  procedures  and  results  im 
some  unsolved  problems. 

Chapter  I  considers  the  problem  of  selecting  a  subset,  containin': 
all  populations  better  than  a  control  under  an  ordering  prior.  Here, 
by  an  ordering  prior  we  mean  that  there  exists  a  known  simple  m 
partial  order  relationship  among  the  unknown  parameters  of  the  treat 
metns  (excluding  the  control).  Three  new  selection  procedure'  are 
proposed  andstejied.  These  procedures  do  meet  the  usual  reuniremmi 
that  the  probability  of  a  correct  selection  is  greater  than  or  equal 
to  a  pre- determined  number  P*.  Two  of  the  three  procedures  use  the 
isotonic  regression  over  the  samnle  means  of  the  k  treatments  with 
respect  to  (wrt)  the  given  ordering  prior.  Tables  which  are  necessary 
to  carry  out  t.he  selection  procedures  with  isofnni<  appniac  h  lor  t  ho 
selection  of  unknown  means  of  norma!  nopul  at.  ions  and  gamma  pupulat  inns 
are  given.  Monte  Carlo  comparisons  on  the  performance  ot  sever, 1 1 
procedures  for  the  normal  or  gamma  means  problem  were  i, irrird  oat  m 


several  selected  cases;  these  are  cj i ve n  in  Table  V  and  Table  VI  at  the 

end  of  Chapter  !.  In  each  case  ten  thousand  simulations  were1  performed. 

The  results  of  this  study  seem  to  indicate  that  the  procedures  based  on 

isotonic  estimators  always  have  superior  performance,  expecially,  when 

there  are  more  than  one  bad  populations  (in  comparison  with  the  control). 

★ 

Chapter  11  deals  with  a  new  'Bayes-P  '  approach  about  the  problem 

of  selecting  a  subset  which  contains  the  'best'  of  k  populations.  Here. 

by  best  we  mean  the  (unknown)  population  with  the  largest  unknown  mean. 

* 

Ihe  (non- randomized)  Bayes-P  rule  refers  to  a  rule  with  minimum  risk  in 

the  class  of  (non- randomized)  rules  which  satisfy  the  condition  f.nnt  trie 

★ 

posterior  probability  of  selecting  the  best  is  at  least  equal  to  P  . 

•A 

(■iven  the  priors  of  the  unknown  parameters,  two  ‘Bayes-P  '  subset  selec¬ 
tion  procedures  and  (randomized  and  non-randomi zed ,  respect. ivel  ,  ' 

under  rertain  loss  functions  are  obtained  and  compared  with  the  classi- 

M 

cal  max  1 mum- type  means  procedure  .  The  comparisons  of  the  performance 
of  with  f.|L  and  based  on  Monte  Carlo  studies,  indicate  that  the 
procedure  always  has  higher  'efficiency'  and  smaller  expected  select- 

n 

ml  '.ize  of  the  selected  subset.  Also  ,  '  appears  t.o  he  robust  when  the 
true  distributions  are  not  normal  but  are  some  other  symmetric  dislribj- 
t  ions  such  as,  the  logistic,  the  double  exponential,  Laplace,  and  > he 
gross  error  model  (the  contaminated  distribution). 


CHAPTFR  I 


SELECTION  PROCEDURES  FOR  POPULATIONS 
BETTER  THAN  A  CONTROL  UNDER  ORDERING  PRIOR 

I . 1 .  IntroductJ on 

In  this  chapter,  three  new  selection  procedures  are  given  fin  i  in- 
problem  of  selecting  a  subset  which  contains  all  populations  belter 
than  a  standard  or  control  under  simple  or  partial  ordering  prim. 

Here  by  simple  or  partial  ordering  prior  we  mean  that  there  rr.i 
known  simple  or  partial  order  relationships  (defined  more  spei  itii.uMy 
later  in  Section  1.2)  among  unknown  parameters.  The  procedures  .in¬ 
scribed  do  meet  the  usual  requirement  that  the  probabi  1  i  t.  i  es  o'  a  (di¬ 
rect  selection  are  greater  than  or  equal  to  a  predetermined  number  I’*, 
the  so-called  P*  condition. 

Many  authors  have  considered  the  problem  of  comparing  populations 
with  a  control  under  different  types  of  formulations  (see  Gupta  .uei 
Panchapakesan  (1979)).  Dunnett  (1955)  considered  the  problem  ol  .«  i ■ 
aratinq  those  treatments  which  are  better  than  the  control  Iron  llio.e 
that  are  worse.  Gupta  and  Sobel  (1959),  <Gupta  (1966),  Naif  (Id/".', 
Brostroin  (  1  977)  studied  the  problem  of  selecting  a  subset  lontaming 
all  populations  better  than  the  control.  Lehmann  (1961)  disies  .  d 
similar  problems  with  emphasis  on  the  derivation  of  a  rest c it  t ed  mini 
max  procedure.  Kim  (1979),  Hsiao  ( I n/q )  studied  the*  problem  o! 


A 


f 


select .  imi  populations  close  to  a  control.  In  all  these  papers  it  is 
assumed  that  all  populations  are  independent  and  that  there  is  no  in¬ 
formation  about  the  order  of  unknown  parameters.  However,  in  man',  sit¬ 
uations,  we  may  know  something  about  the  unknown  parameters.  What  we 
know  is  always  not  the  prior  distributions  but  some  partial  or  incom¬ 
plete  prior  information,  such  as  the  simple  or  partial  order  relation¬ 
ship  among  the  unknown  parameters.  This  type  of  information  about  the 
ordering  prior  may  come  from  the  past  experiences;  or  it  may  arise  in 
the  experiments  where,  for  example,  higher  dose  level  of  some  drugs 
always  has  larger  effect  (side-effect)  on  the  patients. 

fn  Section  ).?  definitions  and  notations  used  in  this  chapter  are 
introduced.  In  Section  1.3  we  consider  the  problem  for  location  pa- 
meters.  We  propose  three  types  of  selection  procedures  for  the  cases 
when  the  control  parameter  is  known  or  not  known  (the  scale  parameter 
may  or  may  not  be  assumed  known).  Some  equivalent  forms  of  the  pro¬ 
cedures  are  given,  and  their  properties  are  discussed.  In  Section  1.4 
the  problem  for  scale  parameters  of  the  gamma  distributions  is  consid¬ 
ered  and  three  analogous  selection  procedures  are  propose!.  In  both 
Section  1.3  and  1.4  simple  ordering  priors  are  assumed  and  some  theo¬ 
rems  in  the  theory  of  random  walks  are  used.  In  Section  1.5  a  selec¬ 
tion  procedure  is  given  for  the  problem  of  selecting  all  populations 
hotter  than  the  control  under  partial  ordering  prior.  Section  l.f 
deals  with  the  use  of  Mon te  Carlo  techniques  to  make  comparisons  amotv; 
the  selection  procedures  proposed  in  Section  1.1  and  those  in  Section 
1.4,  respectively. 


1.2.  Notations  and  De  f i n  i  t i on  s 

Suppose  we  have  k  +  1  populations  >y,  j,...  The  popu I.:  i inn 

treatment  -r  is  called  the  control  or  standard  population.  A-.-.- in" 
that  the  random  variables  X..  associ ated  wi th  F  ( - ;  <  .  )  and 

ij  in 

i  =  l,...,k,  is  an  independent  sample  from  ..  Assume  Mint  we  li  re¬ 
ordering  prior  of  first  we  assume  that,  the  oroerjn;  p 

is  the  simple  order,  so  that  without  loss  of  generality,  we  may  r 
that,  Pj  ...  a.  .  In  Section  1.5  we  will  consider  the  partial  or 
inq  prior  case.  Note  that  the  values  of  “.‘s  are  unknown. 

Suppose  our  goal  is  to  find  a  non-tri vial  (small )  subset  whi.n 
tains  all  populations  with  parameter  larger  (smaller)  than  the  mi 
Pq  (known  or  unknown)  with  orobability  not  less  than  a  pivn  v : ! . e 
The  action  space  C.  is  the  class  of  all  subsets  of  set  I. 

An  action  A  is  the  selection  of  some  subset  of  the  k  population-  . 
means  that  r.  is  included  in  the  selected  subset. 

Let  ■’>  (  -q,  ••  | , .  . .  .  Then  the  parameter  space  is  det.nt.-d  li 

k  + 1 

where  :  =  f"  C  K  !  •  A„  ...  •  ■>„  •  ■  is  a  -..ilr.et 

1-2  k  0 

k+ 1 

k  +  1  dimensional  Euclidean  space  !R 

The  sample  space  is  denoted  by  where 


A.  ~  1  x  t 


n ,  -F  . . +n 
R 


k 


Definition  1.2.1.  A  (non-randomizeri)  selection  procedure  (rule' 
is  a  mapping  from  to  u. 


f 


A  population  " j  ( i  1 . k)  is  called  u  pood  population  if 

"j  "q,  and  we  say  a  selection  procedure  '  make  a  correct  selection 
(f'i)  il  the  *. (*!(■(  toil  subset  contains  all  «|ood  populations.  A  sole, 
tion  procedure  1  satisfies  the  P*-condition  if 
PjCSb)  :  P*  for  all  o  € 

t  ha  t  i  s 

inf  P( ( CS i s )  •  P*.  (!.?.!' 

0  c 

l.et  i-  ■  1  iS  j  i  n  f  P  ( C  S  ]  ■' )  -■  P*  •  be  a  collection  of  all  selection 
"C 

procedures  satisfying  the  P*-condi  tion. 

In  the  sequel  we  will  use  the  isotonic  estimators  (see  Barlow, 
Bartholomew,  Bremner  and  Brunk  (1972)).  Hence  we  give  the  following 
in  if  ions  and  theorems. 

Definition  1.2.2.  Let  the  set  *  be  a  finite  set.  A  binary  relation 
"  "  on  J  is  called  a  simple  order  if  it  is 

(1)  reflexive:  x  ••  x  for  xCJ 

(?)  transitive:  x,  y,  z  tj  and  x  •-  y,  y  ■  z  imply  x  ■  z 

(!)  antisymmetric :  x,  ytj  and  x  •  y,  y  ••  x  imply  x  =  v 

(4)  every  two  elements  are  comparable:  x,  yt...  imply  either 

:<  -  y  or  y  x . 

A  partial  order  on  J  is  a  binary  relation  "  on  J,  such  that  if: 
is  (1)  reflexive,  (?)  transitive,  and  (3)  antisymmetric.  Thus  everv 
s  iiiifile  order  is  a  partial  order.  We  use  poset  (j  ,  )  to  denote  the  sc 
that  has  a  partial  order  binary  relation  "•  "  on  it. 


.’fl 


( !  \  q '  ( t )  t  ( t )  v  (j  ( t )  ma  x  ( •’  ( t ) ,  'i  ( t ! ) . 

Hhmi  we  state  the  till  lowin'!: 

Ihenrem  I. The  class  o  of  a!)  isotonic  function',  on  pose 
a  1 1  hi  vex  cone  and  a  lattice. 

I  here  are  some  algorithms ,  such  as  the  "puol -aq.ijr.ent-vio 
algorithm  (sop  page  Id  of  Barlow,  et.  al.  (197.-’))  or  Ayer,  Bi¬ 
ke  id  ami  Si  1  vp man  (195b)  or  the  "up-and-down  blocks"  a  Inca- i  t 
(  I'tM) ,  whic  h  show  how  to  calculate  the  isotonic  regression  u 
nfiier . 

Ihe  tollowinq  max-inin  formulas  were  qiven  by  Ayer  et.  al. 

I heorem  Id'.  ’.  (max-tain  formulas) 

As  Mine  that  we  have  poset  (..  ;•  )  where  „  I .  ^  .  . 

and  that  (unction  -j:  .  -  Ft,  then  the  isotonic  rearession  u* 

wei  iht  W  has  the  tollowinq  fe.r-nias : 

o*(  .)  ■  max  "tin  Av(s,t) 

1  s  i  f-i 

■  ".ax  min  Av  (s  ,i 

$■  i  t  -s 

s  ifin  max  Avis  ,t  '■ 
t  i  s-  i 

li.n  ;-ia>  A v ( s  ,  t 't 
t  i  s-  t 

win  re 

{ 

< ;  ;  J  \  ' 

'  r  » 

o  /  .  \  r'  1> 


Corollary  1.2.1.  (g  +  c)*  g*  +  c, 

(aq)*  *  aq*,  if  a  0. 

Corollary  1.2.2.  [.  (q*)g  +  <,.(n*)]*  ,  ( cj* ) c|*  -*  <,(q*)»  where 

negative  function  and  c,  is  an  arbitrary  ‘‘unction. 

1.3.  Proposed  Selection  Procedures  ‘or  t.iie  Location  Par, imetc 
To  discuss  some  more  general  results,  we  assume  that  t>'>| .«j 
has  an  absolutely  continuous  location-scale  distribution  bun 


x- 


F(x;  ^ 

)  ■■  r( 

l 

) ,  where  0 

r(xi 

1  for  all/,- 

the  values 

of  V 

and  ,i.  are  unknown,  but  their  or«it*rinij, 

K 

L  is 

know 

n .  Cote  that 

i  n 

this  case  we  rep  1,. if 

parameter  space 

by 

all  other 

(riant  i  t  ies  remainin'!  the  s.n 

Let  us 

define 

the 

subspaco  . 

'k-i  :T>  1-i‘l 

i  -  1 _ k 

-1  and 

let 

subspace  ,, 

.  p  ,  and  s ulr  1 1. : 

k 

n  "  :  ■  C 

;  L 

l 

■  0'  ’ 

then  we  have 

o  ..  Note  that  t.he  , 
i-0  1 

;.q  could  be 

known 

or  i 

inknown.  . 

,  i 

"0 

s  unknown  ,  we  as  ,v  ..  t : 

distribution  of  population  ■  n  is  f tv;  '  )  ,in<1  we  tale  in  ti 

observations  Xqj . from  -  and  t  he  sample  space  tune. 


V“ 
•:  x  ur  u 

,.+nk 

!  x  =  (XQV. . . 

•  X.,  . 

Ulip 

a., i , . . . , x ,  ) , .  Us  int  :  h 

/I  kn, 

k 

t  i  on  1  0 , .  . 

.  ,  k  !  of  parametei 

s  p  a  (,  e 

,  we  h.iv«‘ 

inf  P  (CS!  ) 

i .  C 

-  in! 

1  i  1 

:  inf-'  P  (r:r.  •. )  • , 

i 

for  any  selection  procedure  - 

.  C  JS.  1 

'one  e  the  P*-r  oml  i  t  ion  i  \  i 

:  •.  ,i  i.'>i 

c  Pn il  1 1 ■■ 

I  .1  *  I  )l! 

I 

!  i  nn 

,i 1 1  d 

sav  . 

i  >.  ‘lie 

i1  ■ . 

‘it 

on  'rut 

at  fti. 

■pel 1 1  leu  t 
tu 


■i|i  i  I  V .  i  lent 


t.0 


1  -  r  i  r  j-1  J 

we  have 

inf  P_|I(CS!-s1 )  ==  P(Vi  —  d^j  +  1;k),  i  1 . i.  '  I 


') 

Theorem  J_.  3. 1 .  In  c  'se  I,  (s^  known,  common  known  uni  <  nniiinn  .>  I . 
size  n),  if  d^j  +  -|  is  the  solution  of  equation 

P(V,j  -  -  x)  =-  P*  .  1 .  ■ 

where 


1  ■■ 

V.  =  mm  )  Z  and  1.  are  i.i.d. 
1  I  •  r  i  j  =  ]  J 

i  ~  then  ^  satisfies  the  P*-rond  i  t  ion . 


(-h 


Proof.  For  any  i  ,  1  *:  i  •  k , 

inf  P  (CS  1 , 1  -  P(V.  -  -l;  1  .  ,  )  I'* , 

,  .  ,  1  1  1-1*1  :  l 

1-v-  i 

so  j  satisfies  the  P*-condi  f  ion . 

Therefore,  the  problem  o'  find  inn  the  <l(V's  re<lu<  e  ■  m  t  r  . :  i  m«  :  '  I" 

l  :  k 

ri  i r,  t  r  i  hi  1 1  i  i  ms  of  V  j  , .  .  .  ,  and  .  This  is  ai  h  i  o  v  •  •  ■  1  by  u-  i  n<  i  n,  n  thin 


A 


reins  in  the  I  hem  y  of  rarnlnin  walk. 


>ome  Theorems  in  tile  Theory  of  Random  Vi  <1  1 t 


Suppose  Yj,  V.,,...  a  re  independent  random  variables  with  a  cummnr 
distribution  It  not  ronren  (rated  on  a  half-axis,  i.e.  0  •  P  ( Y  ^  O', 

l’( Y |  0)  ■  1.  The  induced  random  walk  is  the  sequence  of  random 

var i ables 


S0  '  Sn  =  Y1  +--'+  V  11  -  ? . 


1  et 


'»  '  r<s,  0 . Vi  '  "•  s„  o: 


(1. 


and 


(s)  =•-  v  :nsn,  o  •  s  i.  :i. 

n-  1 

then  we  have  the  following  theorem  which  was  discovered  bv  Anderson 
(IdStl.  teller  (1971)  gave  an  elegant  short  proof. 


I  hen  rein  I  .  i.  ? . 


1  ng 


1 

i-  -o 


n=  1 


P(S 

n  n 


0). 


theorem  1.  i.T.  (Teller  (1<»7:)) 
I  et 


I  Iter | 


0,.  .  . 


0)  , 


P(s) 


n 


n 

p  s 
'  n 


'  1 . 


hem  ( 


log  ; a ( s  ' 


n-1 


0) . 


I 


By  symmetry,  the  probabilities 

%  -  P<V  ° . Sn  '  !>■  '  "> 

have  the  generating  function  q  given  by 

log  q(s)  =  l  S"  P(Sn  0).  (1.  :■  I.'j 

n=l  1 

Note:  The  above  two  theorems  remain  valid  if  the  signs  and  .if 
replaced  by  ■  and  respectively. 

Now ,  let 

,  r 

U.  -  max  '  j  -  1,2 .  !  I  •  I  '1 

J  1  t  j  r  i  =  1  1 


V.1  -  min  1  7  7!,  j  1,7 .  (!.  M 

1  1  ■  r<  i  r  i  1  1 

where  7!‘s:  are  i.i.d.  with  absolutely  continuous  c.d.f.  <<(■).  W<- 

would  like  to  apply  Theorem  1.7.7  fo  get  the  distribution  of  Ml  and 
VI,  i  =  1,7 . 

Remark  1.3.1.  The  distribution  of  '.V. ,  j  -  l,...,k  lor  some  1  I, 

will  be  used  whenever  our  goal  is  changed  to  selei  L  a  subset  Ion 

ing  no  population  with  parameter  smaller  than  the  control. 

Theorem  1.7.4.  The  generating  function  g(s)  of  Pfll!  <).  i  1,  . 
i  s 


f 


where 

n 

sn  -  >'(z:  -  x),  n  1, 
n  i-1  1 

if  the  distribution  of  =  ZT  -  x  is  not  concentrated  on  a  half-axis. 

Proof.  Since  the  distribution  of  random  variable  Y.  7!  -  x  is  not 

l  i 

r 

concentrated  on  a  half-axis,  and  Y-'s  are  i.i.d.  let  5  -  V  (Z!  -  x), 

1  r  i=l  1 

r  1  , . . .  ,  k .  Then 

■U\  ■  xi  -  <  max  I  S  O'  =  IS,  •  0,...,S.  ■  O'. 

'  1  rj  r  1  J 

By  leller's  Theorem  1.3.3,  we  complete  the  proof. 

Similarly,  ' V j  •  x  >  =  : S.  ■  0,  i  =  1,  ? . j', 

where 

Si  ~  J  (7-  -  x). 

1  r- 1  r 

I henrem  1.1.5.  The  qeneratinq  function  p(s)  of  P(Vj  ■  x 1  is 

V  sj  P( V '.  -  x)  =  exp  V  1  s"  P(S  0)',  (1.3.  UO 

j-  I  J  n= 1  "  n 

il  the  distribution  of  -  / j  -  x  is  not  concentrated  on  a  hall-axis. 

fornllnry  1.1.1.  Both  Theorem  1.3.4  and  Theorem  1.3.5  hold  for  all 
x  mii  h  t  hat  0  f.(x)  1  . 

i! 


Proof.  Let  -  x,  then 

P( Y-j  <  0)  -  G( x ) 

and 

0  G(x)  ■  1, 

hence  Yj  is  not  concentrated  on  a  half-axis. 


Corollary  1.3.2.  Both  Theorem  1.3.4  and  Theorem  1.3.  b  tin  Id  io>  all 
x  whenever  G  =  c.d.f.  of  N(0,1),  or  f.  =  r  which  is  defined  ,ii  Mm 
beginning  of  Section  1.3. 


Proof.  Followed  immediately  by  Corollary  1.3.1. 

Note  that  in  the  case  of  location  parameter  of  normal  pupal. >t  on 
P(U'  -  -  x)  P(V*  >  x). 

Let 

'■Ax)  -  P(S.  -  0),  j  -  1,  ? . 

J  .1  ,1 


a  ( s ) 


nH 


.  n 

J  , 

n  ft' 


we  have 

p(s)  ,»  s'  P(V.‘  x)  exp  (a(s)). 
:  - 1  T 


l  emma  1.1.1.  p  ^ rl  +  ( s ) 

Proof.  Since  p‘  (s  )  r  p 


V  (nl  n(i:(s)  .... 

i-0 


'  •  a'(s),  'he  result  (an  he  proved  !>v 


tion  on  n. 


I 


theorem  1.1.6.  Under  the  assumption  of  Theorem  1.3.5 


n+ 1 


,» (V •  .  x)  ,  .  1-.  1im  d-.-PisJ 

l(Vn+l  (n+H!sVdsn+1 


}'  P(V.“  •  x).\  . .  .  ,  n  ~  0,  1  ,  ? , 

n  +  1  .f;^  .1  n-.i+l 


( 1  .  -.17) 


Where 


P( ;  X)  1,  VX. 


I’rnul.  Uy  lemma  1.3.1,  we  have 


p(v;,+i 


x)  = 


(n+Y)f 


1  im 

s  '■O'* 


p(n*'>(S) 


jL;0  {ml) !  i !  {n-0 !  p  j>fG)  "1.3 


1  y  P<j)(0) 

n+1  .'n  j!  n+l  -  j 


1  11 

,  >  P(  V x)‘  ,  .. 

n '  1  ()  n  +  l-.] 


larly,  we  have 


P  ( U '  Ll  x)  -  •  P(U  •  •  ,  •  x)  P ( 5  .  0). 

n  +  l  n  +  l  .  v  n-i  +  1  '  i 

l  -() 


( 1 .  s .  1 : . ) 


I.  i.-1.  limiting  Di st ri but  ions  nl  U1  and  V1 

n  n 

let  I  (x)  *  P(U'  x)  and  F  (x)  denote  the  limiting  distribution 
n  n  ■ 

lumt  ion  as  n  •  ■  of  IP.  Suppose  the  distribution  of  random  vari  t.bif 
V |  /j  -  x  is  not  c nnren t ru tod  on  a  half  axis,  then  we  have 


n. 


li 

1  -  F  (x)  =  P'S,  -  0}  +  V  p(S,  0, _ S  , 

n'  '  1  '  1  r- 

r-2 

1  -  F  ( x )  =  1  i  m  t  ( s ) , 

s  -r 

and  apply  Andersen- Pel  1 er  Theorem  1.3.2,  we  have 

F  (x)  -  exp  •-  "  ^  P(S  -■  O'  • . 

rLl  r  r 


Simi larly , 


G  (x)  =  exp  <-  '  P ( S  0) 

r=  1  r  >• 


where 


Let. 


g  ( x )  -  p ( v ;  x). 

G  (-  d|! ; )  -  P*. 


If  Z.,  i  -  l,...,k,  are  independent,  identically  d  i  •.  t  r  jf>ut  fd 
N  ( 0 , 1 ) ,  then  we  can  use  the  roan  rente  formula  of  Theorem  I.  :.n  in 
solve  the  equations  P(V-  -■  -  P*,  i  l,...,l.  H<  m .  in 

Case  I ,  '  •  ( x )  -  : ( -  n j ) . 

Remark  1.3.2.  From  formula  (1.3.  }  we  know  that  d|'j(|.,  ( i  i  . .  .  .1 

does  not.  depend  on  k.  And  w  •»  have  cl!  '  I  .  .  dP:.  Ihca  v  .  i  I  n  >  In 

k  -  H  1  :  k  1  :  i 

k  -  1  (1)  6,  10,  ■  and  P*  -  .99,  ,97r<,  .0',,  .92',,  .9, 

.7,  .63  are  tabulated  in  Table  I. 


Hence  we  have 


!J  r  X .  ,  -•  -  c  . 

j.i  -':k 


U  X.  , 

j=l  J:k 


therefore  the  lemma  is  proved. 


t  .  ' 

J 


Definition  l.j.2.  We  define  a  selection  procedure  1  by  rfp  i  : 
the  inequality  in  the  ith  step  of  procedure  •  ^  hy  the  inei|u.ili 

Xi:k  '  0  -  di:k  J  ’  1  1 . k 

v  n 

where  are  the  smallest  values  such  that  ;  ■.  i!  i 

the  P*-condition. 


Theorem 


d  ^  ^  -  a '  i  =  i  lc 

Ji  : k  ai  : k  ’  1  ! 


The  selection  procedure  ^  and  -j  are  i 


den  I  i  i  a 


Proof,  for  any  i,  1  i  k,  by  T  hoor^'ii  1  ,  ’•!  1 


P*  -  in l  p  ( CS  '  '•  •  =  P'7,, 

t  < 


._,(!)  v 

k  -  i  -t  1  :  k  ‘  !■  -  i  1  !  :  k  ’ 


On  fho  other  h.m.t.  u;inq  the  same  arquments  as  *.e<.t  ion  I.  : 


have 


P*  inf  p  ( i  V  1  '  vi 7 

L  .  r  ' ‘  k  -  i  +  I  :  k 

i 


!  1 ' 

Hence  we  have  d:  /  -  q:  .  ,  i 
l  :  k  ;  :  k 


k-i  •  I  :i 


since  a 


'  1  :  k  i  :  k 


>:.  ,  ,  the  ‘  i  rs  ‘  sm.p  of  and  !  ,in>  identji 


!  i  ‘ 


. .  Wl 


l 


2  , .  .  .  ,k  ,  the  event. 


>v  I  I.  l.P  .mil  Lenina  l.  i.  '.  Hon c t?  selection  proi  eda t'oc. 

in1  i ileii t  j 0,1 1 . 

.  .f>.  Si  hup  Other  Proposed  Selection  Procedures 

In  i  u.i'  I,  we  proposed  some  other  selection  procedures : 

'el  in  i  t  i  i  >ii  1,1.1,  We  define  a  selection  procedure  •  ,  by 

select  •  if  and  only  if  X.  ,  ...  -  ■:  i  -  I . . . 

i  l  :  k  u 

»  n 

’here  d  is  tiie  smallest  value  such  that  • satisfies  the  P*-r.; 


li',nrei!i  l.d.l’.  Under  assumptions  of  Case  i,  md  selection  p>-. 
...  it  we  select  population  .  ,  then  we  will  select  popala^ioi 
or  a  1  1  i  i  . 

cm  1 1  .  S  i  in  e  X  .  X  .  ,  lor  .ill  j  i  . 

i  :  V.  .1  :  k 

v>ilii.it  inn  t  In'  7a  I  ue  !  o' 

lor  any  i  ,  I  i  k.  we  nave 


inf  P  (cs 


k  -  i  *  1  :  k 


I 


r 


by  the  same  argument  tor  selection  procedure  |  and  tier" 


V.  -  Z,  ■ 

i  1  :  i 


Z.  . 


1  r  i  r  1  ' 

We  need  the  constant  <!  such  that  P ' V ;  -  d)  P*  held*,  for  i! 

1  -  i  •  k.  By  Lemma  1.  3.?  we  have  J  -  t )  V .  lienee  we  have  t  h. 

J  I  :  r 

ing  theorem. 


Theorem  1.3. a.  Selection  prr-crduro 

<0) 


s  ci :  i  *  f  i(",  f  he 


with  d  =  d 


1  :k ' 


Corollary  1.3.3.  If  ana  S^  .ire  the  selected  subsets  iv.< 
with  selection  procedures  ,  and  .10,  respectively,  then  1  , 


Proof.  Proof  follows  from  Lemma  1.3.?. 


i  i .!  r  c 


Defjnjtjon  3.4^.  The  procedure  ^  is  defined  as 
Ste.o  1.  Select  tt.,  i  ■  1  and  stop,  if 

d  o-V;/ 

otherwise  reject  - ,  and  go  to  step  2. 
Step  ?..  Select  - .  ,  i  ?  and  stop,  if 


y  a 

?  0  7 


otherwise  reject  ~  and  cio  to  step  3. 


olh 


Step  k-1.  Select  . ,  i  k  -  1  and  stop,  if 


Vi  "n  '  t!k-i  ■  ’ 
*  n 


0 


i>t  herw i so  reject 


k-1 


and  <iu  to  step  k . 


n  op  k .  le<  t  ^  and  stop,  i ; 

,vk  ’  0  "  dk  •  ’ 

v  n 

otherwise  reject  ,  . 

k 

Hero  X.  max  •  and  *j  1  s  are  the  smallest  values  such  tnu* 

.1  1,1  i 

,  satislios  the  P*-c:ondi tion . 


I  va  1  imI  i  on  dI  (1 .  ’  s 

i 


I  or  any  i ,  1  i  k  , 

k-i  +  1 

inf  P  (CS1'-,)  =  inf  P  (  '  1  X  •  rn  -  d  . 

aC  .  !1  3  ;,t  .  !  ,M  J  .n 

i  i 


h'Vl.1  -  “k-M  ~> 


'■('k-ai  ■  -  'Via' 


p*,  z.  r(-;  n,  i). 

i 


I iiis  implies  d,  .  ,  -  d  for  all  i,  and 

k- 1  *  1 

d  -  r](l-P*), 

■  P  (P*),  if  T  is  symmetric 
■_1  (r**) .  if  x.  x(.  r  ?/n). 

‘.imilar  to  the  selection  nroredure  ^  ,  we  have  the  followin']  'he 

Iheorem  l.i.lo.  Selection  procedure  satisfies  the  P*-condition  wi 
d.  r'(M’*). 


wa-,  ciivi'M  ! ’ v  Gupta  .in-1  Sobol  (l')hr'. 


I  he  fill  1 1 >w  i  1 1  ■  |  prm  Oil'll  o 
without  assumimi  any  ordpiin.;  prior: 

Ho!  mil  ion  I.  <.(>.  (ho  selection  procedure  ■  „  is  do*  in-- i  ,j  .  ‘olio  w  ; 

„ :  holm  t  if  ,iml  only  if  X.  .  -  d  i  1.  ,  t 

'I  l  l  0 

.  n 

whom  d  is  tin*  siiiti  1 1  t  constant  sued  Unit  ‘  -  satisfies  u,o  ( :  *  -  <;  oi ,  .1  i 
It  was  shown  tha  t  the  value  d  is  determined  hv  'no  on  i.itioi 

i 

I  ( -  d )  !  -  d‘k 

in' 

l  td)  -  P*'s  i  f  i  is  s yi'imo f  r i c . 


I .  A  Dual  Problem 

Wo  start  with  the  same  assumptions  as  in  Section  | .  ?. 1  Case  1.  b 
i  haii'io  our  i|oal  to  select  a  subset  which  contains  no  had  popii  1  a  t  i 
1  bo  del  in  i  t.ion  of  a  correct  selection  (TS )  will  now  he  <  ha>  f.-e  t  to  so 
a  subset  that  contains  no  had  populations  and  the  -  c.on  1  i  f  ion  will  t 
dot  mod  based  on  this  now  dotiuition  of  i  orrei  t  selection  (t'.S). 

In  location  pa  rami' ter  case,  this  problem  is  a  dial  probles  o‘  Ha 
ori  i  ilia  1  prohlo'',  namely,  1  soloct  a  subset  whiih  contains  all  i|ood 
populations  und"r  or-tor  i  nu  prior  ar  s,,"  p*  i  on  " , 


lino  i'u •  t  hod  to 
■ '  ill  .tiU'ns 
'"i  -o  i  |'i  '  i  e  i  a 


solve  'hi'  iimlne-  is  tliat,  first,  c  haiuie  t  he 
ant  •  l'<  ■  i  .'hies  to  opposite  j  f]n  ;  then  use  a  pi  ...  < 
'  si''"  '  ho  '  ■  .',,.ii''s  all  "now  unod"  populat  tons  . 


where  the  "new  pocH"  popul.it  ions  are  the  "uM  ti.i.t"  ponul.it  ions  <•••!( 

chanqinq  siqns;  finally,  reject  the  selected  subset  .m  l  ke..p  . . 

maindors  as  the  desired  selected  subset,  let  i  l,  1  .  !,. 

! 

the  above  procedure  which  corresponds  to  ■  . ,  i  1,  m- 

t  i  ve 1 y . 

Theorem  1.3.1.;.  The  selection  procedure  . ,  i  i  i  •  , 

.  i 

the  P*-''ondi  t  ion  in  which  the  c- erect  sele.  cion  (rs)  i-kmus  *  mu*  it 
selects  a  subset  which  contains  no  bad  population. 


Proof.  Given  P*  and  obse>’va  t  iors ,  ‘or  any  srleitiiui  ie  ...  ede, 
i  -  1,  P,  ■>.  T ,  after  chanqinq  the  siqns  ot  all  .r.s.n  i.i*«  1 
the  probability  that  the  selected  subset  on  tains  all  new 
populations  is  not  less  than  P*.  If  we  » eject  »he  s<  'a:  .1 

S,  then  the  complement  subset  S'  or  c<  contains  ,mv  'new  .e  •  l  e 
ulations  with  probability  less  than  1-1’*,  Put  tin  "n-w  '  ini' 
latinns  are  the  orininall/  bad  population1  so  wh,  *  w<  (..,•■  i 
the  subset  s'  contains  any  orici i rial  ly  u.nl  popnlut  inn  wrh  pi  1 
less  than  1-P*,  in  other  words,  subset  Sf  t  ortaiir.  no  pup,! 

with  probability  (treater  than  or  equal  to  r* .  1  im  e  'in  :  • . 


all 

arbi trary 

true 

ronfiqui 

t  i  t;n r. 

,  we  tiavc  o'  ole'  n,| 

Remark  1.3..;. 

It 

is  eas,  t 

'»  fs  (  m  •  ‘ 

'  ha  *  '  he  vu  in*-  d 

i  :  t 

was 

used  by  . | 

in 

t  he  i  t  *i  s 

*“r  ;  '* 

!e‘er  lined  ! .  ,  •  On  1 

P<V 

i  + 1  1 

*  i  : f  '  1  1 

i  t  i  '  1  i  i  1 


I  '  !  Iir  >1  i  ■.  1  r  i  fm  t  i  on  I  i  \  symiiet  r  i  ( 


t  hen 


(I.  ,(.  ,)  "  riP 
i  :  k  1  i  :  k 


i  ■>  \ 

Some  Proposed  Selection  Procedures  'V  ,  i  - 


When  is  Unknown 


I  ,)M 


unknown,  (  omon  ‘  known,  common  Stiiiiplo  i.’e  n. 


Pel  mil  ion  I.  We  define  <)  selection  procedure  :  h.  r 

t  tii •  inei|ii  1 1  i  t  ies 


X.  .  -  dM  -  ,  i  1 _ 1 

I  :  1  d  l  :  k 

»  II 


m  in'll  educe  ,  ( Oof  in  i  t  ion  1.1.1;  with 


Xj  -  (r  Xq  -  d  P  ,  i  1 . k,  respect  i  vr  1 ' 

.  n 


n  (-M 

tli  i  * •  f()  >  X^./n,  dl  ^  ,  i  1 . k  tiv  the  sri.i  1  lest  .  oust. 

i  I 

( 7 ) 

th.it  the  seliutiiin  pro*  edun  satisfies  t  h<  •  i'*- .  midi  t  ion 


■  t'il.m  to  the  !,r  e  i,  wv  n.iv.-  tne  tnl  low  in  ;  theori 


,  tl'  •'  '  I  *  •!'.  I  .  .11.  1  i  II  ,10/  I  , 

■  I  I  it  I  i  ill 


I  I  1 


IS  do' i 


,(Vi  f  -  ” 


It  is  easy  to  see  that  :  .  ,  ,  -  d;'!.  I'he  t u  1  low i m 

k- i  + 1 : k  I  :  i 

qi  ves  us  an  identical  form  o f  the  selection  p^ncedur  e 


Theorem  1.3.15.  The  selection  procedure  J  will  no?  ho  .  h.n.  n  ! 
the  statistics  X.  ,,  i  -=  1 . k,  arc  replaced  by  X.  ..  i  I.  j 

•  l  K  1  !  r 


respecti vely . 

Proof.  The  proof  is  the  same  as  that  in  Case  1  and  hem,'  i<  i 
ted. 

The  values  d^:,  i  -  l,...,k  are  tabulated  in  Fable 
k  =  1  (1)  6,  B,  10,  and  P*  .99,  .  975 .  .95,  .9P5,  n. 

.75,  .70, .65. 

Similar  to  the  Case  I,  we  propose  a  selection  procedure 
fol lows : 

Oef ini tion  1 . 3. 8.  We  define  a  selection  procedure  \  ^  by 

'■  Select  if  and  only  i:  X.  ^  xn  -  d  i  I,.. 

»  n 

/  ?  \ 

where  d  is  the  smallest  value  such  that  ; .  satisfies  the  1“  . 
Then,  similar  to  Theorem  l.0.0  we  have: 

Theorem  1.  T.16.  Under  as sump'linns  of  ;.,i  ><■  ,1,  t  ?i< •  ,elm  ;  i.»- 
j  satisfies  the  P*  - 1  ()<•  I  i  *  1  ,r  wiM  d  , !, 


fle/.t,  wr  defire  a  selockior 


■  m  ra  ,'d'irf ■  '  whii  h  is  m  i 


but  re  *• :  I  ,ic  i'  by  •  th,-  sr  pie  '•i',au  o'  ’  >«  ■ :  n  j  1 ,  r  t  <n 

0  i)  U 


Me  I  m  i  f  i  on  I  . 


I  tic  so  I  <•(  t  i  nn  procedure 


i  '  ,1c  ’  i  III'I !  ti  /  I  C! ,  i  II  I  I 

X.  -it.  in  .  (Pefiiii  tion  l .  M!  bv  X.  -  <t|  .  1  ’  •  • 

i  i)  i  i  1  n  i 

.  n  • ' 

where  it*  .  .  ,<l.'  tire  the  sin, it  lest  v.i  1  nos  such  th.it  satisfies  'he 

!  k 

!'*  ,  , hi.  1  i  t  inn. 


Similar  to  Thenren  1.3.  1C  we  have: 

I  lieorem  I  1.1/.  Itie  '.election  procedure  satis* 11  ■  *- 

with  ,|.  d,  i  l,...,k  where  d  is  determined  bv  the  rip.,,!t  ism 

/  fl  -  I (t-d) jdC(t)  •  P*,  ’  •  ■ 

/  I  (it- 1  )dl'(  t )  -  P*,  if  ;  is  syr’ft  ri  < 

Ami  *  will  not  he  chanped  if  the  statistics  <■  is  f{?i>la.  e  i  bv  ■  . 
I  he  -..n;i|  1 1  *  *  moan  of  nonuiation  ■  for  i  k. 

!hr  followin')  selection  procedure  '  was  proposed  by  Ounta 
'.libel  (I'thfl): 

Definition  1.1.10.  the  selection  procedure  ,  '  i  defined  bv 

'  •  So  1  a  t  i  1  and  onl  v  i  f  X  •  X  -  I  i  1  . .  . 

1  '  1  10 

i 

where  i|  is  determined  by  the  tnllnwinn  equation  ir  1  is  uee'-al  1 1 
i  *'  i  Ini  I  i  on  : 

■  k  /  . 

■  [F  t  /  a  d)^  (uMu  -  P*.  '  1 .  . 

lit  / 


4 


Tor  the  special  case  n.  -  n  M  ),  I  . .  .  A 

1 

■  !kl  r  +  J;f  i'\ 

If  F  is  non;;,'il  .  i  i  s  t,r;  hut.  ion  N(0,l),  the  tables  <>»  1 
isfyinq  the  Equation  f]..i.?5'  *nr  several  ;alnn*  it  I1'* 
Uechhofpr  ( 1  954  )  *  >r  •  I  1 ' ,  ;  n  s  in  ,|  *  ,i  ■'  1 1  i 1  f . ; 

/  \ 

1 .  .5.9.  Sone  Promised  :.(*]•'*(  t  '•»  :Vo--e  lares  ‘ '  \  \ 

l 

When  l.Diiii’ion  Variate  (  is  hnln ws 
Case  III.  ( ^  know' .  I  t/H  m  van  '■ . •  n < *  :  4  inwti ,  n  .  1 . 

In  this  i  ii'e ,  w>  asv.-ie  tin*  •  \  v.  h  i  ■  '•  i  .  th 

N  { 0 . 1 ) . 

Definition  !.•:.!!.  Wn  define  t  he  ■,..!<  .  ‘  i  m  jm.,,,  !■,'.•  ' 

’  In-  i  I  i  •  i  i", 

■Vt  r,  -  -V:'  1M  1  1 . 1 

in  ;■  fit<  •>,!•, ,,■  |  (  ?'•  ‘  ;  r ■  i  t .  i  1 1 ■  I  .  ’ .  I  '■  , 

’  •  .  -  ■ '  i  i  ’  .  .  .  .  I  .  i  ' " . :  1 1  ■ .  t  i  v  i  ■ 

l  :  V  r  •  :  I 

, !  0 

where  4  ’  S  I’V  the  Si::.l  ,  'e',»  v.t  '  •  II  ,*  : 

I 

r’*- 1  i  in 1 1  i  t  i  on  ,  V  denote  •  t  •  ••  ■  1  >  ,»  .  •  '  !  m  . 

k  'll  -  1  !  ,  that  is 


i 


Null-  l  lt.it 


h.r,  |  hi-  i  hi  -‘,<|ii.irc>  ili-.t  rihut  iuu  with 


t  reeiloii,. 

Ily  ir,  ini)  similar  arciuments  as  in  Case  i,  we  have: 


f 


Remark  1.3.4.  The  values  of  cl'  ,  ,  .  i 

k  -  n  i  :  k 


,  V  - )  i  ■  J :  *  ■  > :  1 1 


heiu  e  .if  ^  t  ,  .  f  ii ;  ' : . 

k-i  H  :k  1  :  i 


B  :  usiri'i  Rabimwit/  and  Weiss  table  { 1 9f»9 )  ( w i  1 1 . 


k-  i  > ; 
.  ’0,  .  ,’b ,  . 


we  have  evaluated  and  tabulated  the  values  of  t, 
in  Table  III,  for  k  ?  (1'  6. 

.80,  and  .  7i>  with  common  savole  si.-c-  n  s,  k  it,: 

for  k  (>  and  n  81,  wo  ■  o  i.  ■  ! ,  1  .  as  an  a ;  I"  n 


i  !  .  . 


I lif*  proof  of  t.  ho  following  theorem  uses  the  same  n  r- •  > j  .>•■.< 


r,  t h.t t  in  C.iso  I,  honco  it  is  omitto'l. 

theorem  1.3.  19.  the  equation  which  determines  tho  tor.stant  3  of 
iHcttion  procedure  is 

t  ■■(yd)ci  (y)dy  -  P*.  -i. 

0 

(iupt.it  <in<1  Sobol  (199H)  gave  a  selection  procedure  •  ij  ’  in  ft; 
, i ■ . i • .  If  is  as  follows: 

[  :  Select  ■  if  and  only  if  X.  ...  -  .1  '  i  1,..., 

1  i  i  0  ■ 

.  n 

l 

im!  the  equation  which  determines  d  is 

/  : k ( yd  in  f y)d/  -  P*.  (1. 

0 


(4 ' 

I.  .10.  1  nine  Proposed  Selection  Procedures  •:  ,  i  ■  1.  3,  -1 

When  iloth  Control  ...  and  Common  Variance  •'  are  Unknown 

ire  IV.  unknown,  common  variance  unknown  and  <  onisrr:  samp 
size  n. 

Wo  av anno  that  in  this  ,  ase  distribution  f  is  tho  c.il.t.  f; ; f i 


n 


( 


f  •  ( V . 


,(4' 

“i.  -  |  a  !  :i 


f  ■  ■!.  /  ; 


The  constant  .1  of  procedure  V'  is 

,  ..  i!4) 

'  '  v1:k- 

rho  cons* ants  i  of  nrocedures  •  t<n.J  !  arc  iletei !■ 


■  f  :r(.j  -  tdV  ' tj ^ .1 . •  ( ♦ 
0  - 


with  r  -  i  .?n ■  i  k,  respectively,  and  their  values  fur  sole*  i<-  .  v.i!  i 
of  P*.  k  am)  are  niven  in  Gupta  a"1  Sr.be  I  ( 1  ')[■  •'  1  and  •  .nu'-t  \  ; 

1..1.11.  Prupi  rtiu  (  [  tp,.  Se  1  pc  *  i  or  Prut  e lur.  s 

Under  simple  ordorimi  prior,  it  ;s  natural  to  i , ii  <■  in.r  m  > 
selection  procedure  is  O'-der-preservin.j  ar  define!  below: 

Definition  1 .  >.14.  A  select  ion  r>r<v  <■  turn  is  m  ler-nre' i.v  in  ,  j  • 

selects  .  w  i  t  r  parameter  .,  arx  i '  then  i  •  a!  .>■  i  ! .  .  • 

i  :  i 

Procedure  is  weak  order-preservin':  ,,r  ;t 

P(  is  selected'  )  r'  .  is  sr'rrted  )  whenem* 

l  l  '  i 


It  is  en  :  •/  to  see  that  any  f  ?  »~i  If  ■ »  -  pi  esei  v  i  n-;  sole  t  i . p  r  k  .  ,  i  ,  i 

is  weak  order-preservin':,  bm  ‘dr  <  .  .  v<  r  a  i  >...t  ‘mi  . 

1 1  ) 

flow .  let  ■  :'>  -  i  1  t 

1  l 


Theorem  l.’.r’O.  The  selection  proi  e  .  an 


i  i  ' 


in' 


order- ;  re  $  o  r  V  1  no  and  prnc  (‘ilijre 


1 


1  )M.  :  t  I  I'll  ■  .  f  1  ,  l 


On  the  other  hand. 


or  iirocodiire 


*  •  ■  ■  i. 

i 

"o 

i  + 1  ' 

,  i  f 

c 

i 

■  "k 

0  ' 

i  f  i 

-  k. 

-  -r-  e 

i 

'  0  ’ 

"l f’ 

i  f  i 

-  0. 

1.4.1.  Proposed  Selection  Procedures  i  f>,  /, 
Case  I.  Control  ^  known  and  common  sample  si/e  n. 

Definition  1.4.1.  The  selection  procedure  is  -let 
Step  1.  Select  i  k  and  stop,  i f 

Y  r 

“k : k  '  V  :  k  0’ 

otherwise  reject  ^  and  no  to  step 
Step  ?.  Select  i  •  k  -  1  and  ston,  it 

Xk-1  :k  '  Ck-1  -k  '  O' 

otherwise  reject  ^  ^  and  no  to  step  ■. 

Step  k-1.  Select  i  ,  .n'.d  s’ op,  it 

v  .  , 

,  :t  O' 

otherwise  reject  „  -trip  <|e  to  step  l. 

Step  k.  Select  ^  and  stop.  i‘ 

''  1  :  i  '  ’  :  k  O’ 


otherwise  re.iei 


I  '  ,  I  |  ,  .  .  .  ,  r  .!!•••  !  Pi-  .,1  I  1  !"■  ' 


l  n‘i  m  ,•  luri-  ..I  f  i  .  ’  it"  t  >|.«  01,  I  I  '  MU, . 


;  hi'iit', '.i  1.1.1.  \/<-  h.ivt'  i.o  ii  , ,n  ','i  ,t)U>  '  i,’i 


i  |  ,  ‘inA  f  ho  i  oust  ,ir.  • 


1  i it •  ■  i 


f  .  ) 

'  i  i  M 


Y  V 

•t*  (  '  tl  ' 


i  v .  ,i  re  i .  i .  A.  with  ,leir,  i  t  y 


•  i (  •  ;  ,  ,  ,  .hi, ;  ;  .  1 .  f  .  , ( 


i  h,,'i  Mir  ;,rn( .  crliiff;  ‘..it  isfio’,  Mu-  '•'*  - ,  on  .1  i  *  i  oh  . 


•l  ,1,  ,  1  .  Mil'  ‘III  7  I 


f  1  ho  m Mi, iw"  •  i  ;i ■ 


fli'ii'  ill-  i  ‘|(io  1  |io|iii  1  ,i  t  i  on ,  Mirn.  .ti'-'i'r  Mio  j ,  *-f  >: 


inf  (rt 


k  Y  +  ...-•  Y 

inf  P  (  LI  (  'iia  x  mi  n  s  -  s  .  , ,  (  .  .  . 

,  .  .  .  ,  t-sH  i  : I 


=  LI  (  max  min 

M  1  < s -  i  j-t-i 


Ys  "0  +---+  Yt  n 


t-s  +  1 


- 1  i 


/  ♦...*•  Y  - 

t  s  i 

,  max  -  v  ■  *  c  •  .  1 

,  .  l -s  +  i  -  i  :  k • 

Is-  1 


-  Pfl) .  c.  , 

‘  l  i  :  k ' 


where  Y.'s  are  i.i.d  with  density  of-;  ), 

i  J  ti 


.  i 


and 


U. 

l 


Y 

max 
1-  s-'i 


a  .  +  Y  . 

l 

i  ->,  *  1 


Corollary  1.4.1.  c  i  -  1  ....  ,k . 

-  •  J  i  :k  i  :  i 

M 

for  any  *  0,  let.  -  f’v.  -  / 1 .  n  I,  ' 

n  .  j  i  M 

0  •  P  i  <  ■  -  /  0  )  1 1 1  x  ,  .  ,  )  I  ,  '  In  ■  1 1  i  .  t  r  i  I  ‘.i  I  i  i  in  i  < !  Y 

t  n  I 

(  oni.en  t  ra  t.eM  a  ha  I  )-•;>.!  .  ;  /  "h<-  h'im  i  ].  .i  tin-  innhal'i! 

inn  toon  *  ji  ri  <>•  i  'i'"'ilat.  i  ve  di  S  ?  r  i  foot  '  on  taint  inn1.  I'flJ^  ■> 
-  I  ,  .  is  niven  by 


k 


Mt'iii  t ■  hy  Ihrorei'i  1.  I.D,  wo  i.ave  the  followin'!  i'i'nn  rn  i 
ill  <  U 


''("mi  *> 


r  i  '<  t‘(  ■  ( .i  *  1  )ti ;  :  j.|  j  .  I 

kn  I  0  k  1 


where 


>(U 


0 


When  x  0,  both  si<iec,  of  I  q.iation  (1.4.4)  equal  to  /ei  o 
. i  I  ■ . 1 1  ho !  >(•.  f  or  x  (). 

Note  tint 


I 


P(  .  V  Y  .  •  x)  ■-  !\  ( x  ;  r  , 


r  i-1  1 


rn 


■  u  (  xrri ;  r  ,  I  ) . 

the  values  ck.k(P*,  . ,  n)  which  satisfy  f  quati>>n  '1 
t ,  1 1 1 1 1 1  a  t  ed  in  Table  IV  for  k  -  1  (1)  It),  P*  -  .  )9,  .  % ,  . 1 
T,  4,  fi,  ati'l  n  4,  I. ,  ID,  ib,  ;’pi. 


ei.n'iii  1.4.1.  r  j  ,  k  i  .  (  ,  .  k  f  or  alll  i  k  -  !  . 


Proof.  ftie  .  ore.  f  ari  f />  r  .  _  ^  (l  ),...,()  are  he'eri  p 


i  i " . |  iei  t  i  ye  I  y  . 


i  '  I 


l .  .  0  10  1  i  O',  (  .  .  .  ,  f  or 

lit  i »  ;  ;  k 


1.4. 


I  he  1  ei  *  i  i  r.  arm  eu  m 


VJ 1  1  1  not  tie 


i  lif  I  f '  t  h. m<l  s  i  lie  i the 


■  Iff  l'fCr,  i:!  ill 


f  .  of  '  with 
hfin  c  Mu  vi  loo  *  1  •  fin  bo  o.i  i  1  y  so  I  vod  with  t  ho  ho  Ip 
i  ti  i  -  sipian ■  (1  i  s  t.r  i  hti  f  i  on  . 


1  i  t  .1  t  i 

i  ill 

to  trie  Soli 

n  1 i on  o f 

V-.riam 

P  ()  j 

“  irra  1  [' 

j  i 

i 

/  t  .  i  n. 
1 

,  1 . r 

,|  T'l  tfl 

'  ‘Mil 

lo  par.rio 

»l»  »r" 

ia  1  hoihi  1 

.it  ions  and 

...  ft 

1  _ _ ' 

• ;  i 

1 . k 

•,crMt  ions  on  the  population  with  the  mean  ....  (know.  . 

Mu 1 1  I  ho  order  is  known. 

1  k 

In  the  appl  icnt  ion  of  selection  procedure  ^  or  w.n.M 
'  ii  do  is  to  evaluate  the  isotonic  repression  of  S',  whiv.h  is 

M 

v.ii  iriiKO  of  population  i  -  1 . k  and  denote  t  by  St 

i  hen  dim  My  apply  •  or  The  constant  we  neel  is  let  ft 

o 

I  .pint  ions  (I.A.?)  and  (1.1. d;  where  we  replace  bv  <»,  the  i 

»  n 

Mull  ’  n  /  .  lias  .'  distribution  will  n  degrees  of  ?  reodo1  < 
i  i  n 

i .  '  ■  !  .  has  the  r . d . I . 

i  i 

r>(?nt;  n,  ?)  ‘  G(  t  ;  n,  ^ '  , 

hr  in  f 

|  *' 

P '  V.  ft  (i(  t  rn ;  rn,  ]). 
iM  1 

I  hf*  ipplirjtion  -v  .  i‘.  i!?'iiar  !n  t  M«i  f  <V  u  . '  pp  Mr1 

(  I' W  hri  t  wp  nor*'!  t»  -Ip  in  *n  I'crl.n  7  .  in  v.  :  :j  i 


T  II*- 

.P  i:V 

:  .*»  ( *  r 
thp  1 
■i 

'ml  HP 

'OiiSCV' 


Remark  1.4.1.  Q  'G..pta  and  'lobe!  ( 1%:  )  5  ■  to* -  .  not  l<- 
inq  prior  and  the  sample  sizes  fnr  each  pojtuiat  i •  >f •  ••e<- 


It  the  ■■leans  ....  1  1 . t  are  un  kriov.'ii  in  1  » ,*:<s  ■ 

l 


n  -  1,  let  S ;  (X.  ■  -  v. .  ■  /n- 1  and  :>so  n  -  i  in 

l 


.  A  . 

1  n  1 


(1.4.4),  (1,4.6)  and  M. which  de t er-.i ned  the  ..or 

c  and  c  1  for  ,  jr.i  resnec  *  ,:vo1v. 

c  11 


1.4.?.  Selection  Procedure  1  ,1  6,  /,  o 

i 

Case  II.  p  unknown. 

The  assumptions  ire  the  s-vo  as  in  '  .isi-  !  ept  t  n  • 1  n 
tioris,  viz.  ,  xm  , .  .  .  ,  k„  are  taker  on 

11 1  Urip  u 


i 

Tor  selection  procedure  g  the  inequalities 
cedure  ar-.d  rorrespondinq  to  X..,  c'-k  o  ■  ^ 


■  if  tin 


.  , .  ,  r  |  •  I 


c ) ( i  1 . k ) .  respecti  \ 


i  :  k  '  i  •  k  ‘ ' ) 


I  ti. 


..(2)  ,• 


A 


i  .  obtained  as  before  v 


riven 


f  p;i!,  i  .  ‘  i'  :  f  !  :t 


i  }  | ,  , f  • 


w hem  pf m  ■  is  the  tv  i- 

0  .  d .  ‘  .  of  /  ■  I '  POliii  1 .1  t  i  S’l 


If  population  p  has  q.rr-a  di  s  f  r  ■  b:i  t  i  on  with  dr 


1  in  known  onloown  . 


lor  so  1  r (  r  i  ■  in  ;  ■  >'i  ii  ed1 1  >‘e 


(I  O'  n 


ail 

)  . 

I 


1  n  n  I  , :  I  I  .  I  ■ 


I 


'  I  , 


k  1 II 

. leprii.lent  p. ipu  1  a t  inns,  p.u  n  imp 

u  1  tit.  i 

■  .1.  W  i  t  •' 

•  list  r  i 

!)ii  I.  i  f  > r i  \  ( 

will) 

i  iimiiicii  known  var  i  .nu  i  ■ 

.mil  ( 

i  >!,ii- .on 

sai  a.  1  .• 

.  i  .*(*  1; 

/  •. •. .ini- 

1  In  ■  i 

hip. in  o  1  he  1  on  ti  n  1  i  s 

know?" 

;  without.  In 

' i  !  : 

.  i  ]  i  t  v 

assume  tJi.it  0  .'Hid  -/.n 

i  . 

In  tfie  s  iniii  1  <i t  i on  ,  wt*  use 

•tub  1  IT 

and 

nink i “ 1 

s  LVr'-:‘ 

siiv.iO*--  Vo i' 

I'.M  t 

,n|p,  I’urdiiP  University  Co 

i'UU.  r  ir. 

a  ten 

ter.  fi. 

'ip'.t**-  * 

‘  P  r  e  ■!'  )' 

1 1 1 iiul )f •  i .  fiir  p.u  li  k ,  wo  pene 

ra  ted 

one  r 

andom  n 

.r:\wr  , 

vari  i ! H  p 

imi  h 

l it i[ id  1  ,i t  i on ,  then  app  1  i e  ; 

each 

soles 

t  ion  ra¬ 

or  •  »= t n  *  ► 

see  "'a ;  s 

.  i  n  •  1 

repeated  it.  ten  thousand 

t.  imes : 

we  u 

sed  t  hi¬ 

v'(>  I  ill 

VP  riup.p 

.r.  in  approx imat i on  ot  the  tv 

act  values 

nt  the 

tl  r>  •  •  >;.  1  •! 

ted 

i  har 

at  t  er  i s  t  if  s  f  or-  each  pro.: 

edure. 

In 

Table  V 

we  iM- 

Hr.  Hi!' 

Mill  .1 

1  inns : 

.  ■ ,  -  is  Lne  u.ira 

1  k  i 

.meter 

■  of  | i 

'  i  1  a  v  i 

' 

i"s  I'HO 

PI  P(i  nrrect.  1  v  reject' in 

d  all 

bad  pupal  at  i . hi •- ' 

IT.  I'(  correct  c  lassifka 

t  i  an  n 

. !'  a  1  1 

popul  a 

t  1  dll  ' 

wherr 1  I  tn ■  1  nrrec t  i  1  ass  i  f 

i  '-.l  f  i  i 

r  -ea 

ns  that 

W"  m.; 

*  . ; 

pupil  1  a t  i ons  and  re  jei  1  a  i 

1  ev: 

•■■Ipu  i 

a  tie”  . 

:  1  i  >. pm  t  ed  number  (  ,  i 

e  )  r,  • 

hi»  «  •  UOjMl  Lit  i 

cats  .  i "i 

!  .1  in..  1  '  a 

Hi”  selected  subset. 

‘  ;  -  r/  ' 

i  0 

i 

.  iff* 

i j  ' 

:  1  <pei  t  i'd  S  i  /e  u  *  t . 

, f‘  j  1* , 

■  r  ' 

jtlV-  •  . 

.  Ad  L  l 


l.ition  result. s 

for  the  co" 

t>u  ru  t  'i  ve 

IJt-rtorHiiiK. e  ot  ; 

e  1 "  .  ■ 

•(lures  for  the 

<!<rii']ti  niedii 

pro  1)1  on 

(notiitien  exp  Ini 

K'O  in 
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,  I ,k  ; .  i  action  A  <  o  i  ■-  i  no  st  lection  i 
populations.  i  (  A  means  that  ^  is  included  ii 
fiction  At  i.  is  cal  loil  a  correct  selection  (At) 
is  included  in  the  selected  subset  A. 

Definition  / .  \ .  / .  A  pleasurable  function  del 
select  ion  procedure  provided  that  for  each  xt  ; 
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Lemma  2.2 A. 
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Proof.  follows  from  the  definition 


*  l 

:neore::i  P.<f.l.  Gwen  a  numher  I’  •  2  ■  I),  the  prim 

t unction  L j ,  the  selection  procedure  . ^  is  ,i  non-isindm 
rule. 


Proof.  It  is  sufficient  to  show  that  the  selection  pm  ,-.w 

* 

the  smallest  posterior  risk  in  the  class  ; .  ■  iv.  - 

A 

tion  X  •  x.  Let  the  posterior  risk  of  c  I’  )  hr 

N  K 
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if  *jjR(x)  -  (0 . 0,1,...,]). 

s  terms 


EFF(i|.'  ) x ) 


i  =  k.s+l  ^NRli  )(-)P[i](-)  +  cp[k-s]( 


X  VNR(i ) 


(x)  +  c 


X  #NR(i)(-x)p[i]W 

k  b 

X  ’NR(i)(x) 


=  EFF(^Jx) 


The  last  inequality  is  obtained  by 


X  C(i)(-x>[i](x)  =  jk_s+1  4(i)(-x)p[i]( 


SLs+l  ,"NR(i)(“))  Ptk'S](x>' 


theorem  2.4.4.  The  randomized  selection  procedure  v  is  the  PME  proce¬ 


dure  in  t,  (i,P  )  =  l  for  given  r,  P  . 


Proof.  It  suffices  to  show  that,  given  r,  P  ,  X  =  x, 

EFT  UB|x)  ■  E F F ( '  |x),  \  (  i  1  . 

Suppose  --  (0, . . .  ,  ,1  , . .  .  ,1 )  0  <  v  <  1,  1  <  s  k  -  1. 

s  tenns 


By  theorem  2.3.1  there  exists  c  ■  0  such  that 


k  k 

r  u  \  _  v 


T,  * c  '  T,  *(o(x) 


<)/ 


If  0  <  c  <  1,  then 


EFF(^  |x)  = 


k 

l  * 

i  =  l 


(i)(x)p[i](x) 


I  'F(-i)fxJ 

i  =  l  '  ; 

k  B 

*(i)(x)u[i](x)  +  cp[k-s](x) 


i  =  l 


(x)  +  c 


h  v(i )(x)p[i](x) 


k  R 

Jl  ^ 

B,. 


=  EFF( ^  J x)  . 

If  1  <  c  =  v‘  +t+(l-v),t>0  integer,  0  -  V  -  1  then 
k  k 

f,  ---ci !{*,p[i](x>  *  "'"ik-.-mf 


+  P[k-s-t](x)  +  •••  +  (1-  )p[k- 


<'(i)(x)p[i](x)  +  (p[k-s](x) 
hence  by  the  same  argument  as  above  we  have 

EFF ( v*  ’x)  -  L ‘  F ( ;<B | x )  . 

Since  x  is  arbitrary,  the  result  holds  for  all  x. 

2.5.  Appl icat ions  to  Normal  Model 

Suppose  we  have  k  populations  population  n.  has  distnbu 

2 

tion  N(u^,o.),  where  o.’s  are  known  and 


e-’s  are  unknown.  Assume  that  we 


t 


have  independent  observations  ,  i  =  Let 

n.  1  ^ 

I 

X.  »  X.  .  and  let  X  =  (X,,...  ,X. ). 

i  n  -  ij  -  v  1  k' 

Suppose  we  are  interested  in  selecting  a  subset  containing  the  best 

★ 

(the  population  having  the  largest  mean)  under  the  posterior-P  condition. 

★ 

wrl  some  prior  i  =  i(u).  Then  to  find  a  Bayes-P  selection  procedure  is 
equivalent,  in  some  sense,  to  finding  p^x),  which  is  the  posterior  prob¬ 
ability  of  the  event  ( u ■  is  the  best!  ,  given  observations  X  -  x,  wrt  a 
given  prior  i,  for  all  i  =  l,...,k. 


Case  1.  Assume  that  we  have  a  common  sample  size  n  and  a  common  known 
variance  ■  ^ . 

la.  Suppose  we  have  no  prior  information  about  the  unknown  parameters, 

and  use  the  "non- i nformati ve"  (Box  and  Tiao  (1973))  or  "locally  uniform" 

prior  p(u.)«.c  for  each  population. 

The  posterior  density  function  g-  of  ,i  • ,  given  x  is  the  normal  den- 

21 

sity  with  mean  x.  and  variance  '  /n,  i.e., 

2 

•  n(,,  -x.) 

‘h(";  !*)  "  - -  exp  (-  -  o-  -  )  . 

1  .2-,  •  ?  • 


P(i](x)  P(i-(i)  -  -;k],x  =  *' 


i  ^  ;(t  ♦  (X[i]  -  *Li])d:(t) 


i  1 , . . . ,k. 


Men'  r,  the  quantity  corresponding  to  the  i  largest  observation 


i 


lb.  If  u-'s  are  independent  and  have  the  identical  prior  distribution 
2  2 

N(t>o»°0)  anc*  Xilpi  N(p.,Oj/n),  then  it  is  well  known  that  the  posterior 
density  funciton  g^  of  p.,  given  X  -  x  is 

9 

9i(ui  | x)  -  N(t  ,t,  )  with  SIP  property 


where 


Hence 


%  =  <V’o  +  nai2xi) 


P[i](x)  =  *(t  +  tn,Ji2<x[i]  - 

The  last  expression  for  p^.^x)  is  the  same  as  that  for  the  non- informa¬ 
tive  prior  whenever  Og  -+  ■*>. 

Since  p^j(x)  =  Pj-^j(x  +  b)  and  since  the  normal  distribution  has 
the  strictly  SIP,  it  follows  that  and  are  "just"  a.e.  and  trans¬ 
lation-invariant  in  both  case  la  and  lb. 


Case  II-  Variance  o^'s  are  known  but  j i ' s  and  n. 's  are  not  all  equal. 
1 1 a .  Using  the  non-informati ve  prior  p(p.)  c,  rl . k,  we  have 

ix>  V  i  \  Y  -i  —  V 

P(i)(x)  *  /  a  ■!>( t  a  _W  .tjJ  )d*(t) 

(  ’  -  jJ'l  V(J)  (j) 


a  /  .  v 

where  V/.N  =  -  - 

<’>  "(i) 


i  -  l,...,k.  p^.  j ,  and  n^.^  are  correspondin'!  to 


and  we  have  the  following  theorem. 
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theorem  2.5.1.  p^.^(x)  is  non-decreasing  in  i,  i.e. 


p(i)(x)  - 


>[i] 


(x). 


Remark  2.5.1.  From  the  above  formula  of  p^(x),  it  is  easy  to  see,  in¬ 
creasing  the  sample  size  of  the  non-best  populations  will  increase  the 
probability  that  the  best  population  to  be  selected,  however,  before  do¬ 
ing  this,  we  don't  know  which  one  is  the  best  one. 

B  B 

In  this  case  ij;  and  are  "just"  a.e.  and  i  anslation-invariant. 


Case  III.  Assume  that  priors  are  independent  but  not  identical  normal 

2 

distributions,  namely,  -  N(e^,Og^),  where  e^'s  are  not  all  equal;  if 

2 

°li 

the  conditional  distribution  of  X.,  given  u-,  is  N(ii^,  —  ),  then  the 

'  i 

posterior  density  of  n.j>  given  X.  =  x-  is  g^ujx.),  which  is  the  prob- 

p 

ability  density  function  of  normal  distritubion  N(f>x  , r. . )  where 


\  =  4(o0i°i  +  nialixi) 


Hence  we  have 


p.  (x)  =  /  n  o[t  ---  +  (o  -  T>  )]d$(t). 
1  -•  j/i  f'j  ’j  xi  xj 

If  ,i(),  o^.  =  and  n.  n,  i  -  l,...,k,  then 

f ..  -  f.  =  (o/  +  i  -  l,...,k 

and 


P1  vx) 


/  >Lt  + 

.17  i 


,V"j  + 

\  2  -  V 

0  ’l 


)  ]d<t(t )  . 


10 


Case  IV.  The  General  Normal  Model 

Here  we  consider  a  more  general  prior.  Suppose  we  have  k  popu¬ 
lations,  common  sample  size  n  for  each  population,  and  common  known 
2 

variance  o  ■*  0.  The  observation  can  reduce  to  X  =  (X-j , . . .  ,Xj. )  whore 
n 

X.  =  l  X../n,  by  sufficiency. 

1  i=l  1J 

The  "Normal  Model"  is  defined  as  follows: 

2 

(a)  X| n  '  N(u,  ql),  q  -  ~ 

where  I  is  the  k  x  k  identity  matrix. 

So  the  X's  are  (condi tional ly) independent  when  u  is  given. 

(b)  vi  •  N(0fl  1,  yl  +  tU) 

where  eQ  €  F  ,  Y  >  0,  t  >  -  ■*, 

1  =  (1,...,1)  and  U  =  V  1  . 

Here  y  >  0  and  t  >  -  ^  are  necessary  and  sufficient  for  ,1  t  (I! 
to  be  positive  definite.  This  model  was  chosen  by  Chernoff  and  Ynlinv 
(1977)  (t  •>  0),  Gupta  and  Hsu  (1978)  and  Miescke  (1979). 

By  (a)  and  (b)  we  qet  the  posterior  distribution  of  ,  given 
X  =  x,  and  the  distribution  of  X  as  follows: 

p  1  x  N(‘>,  al  +  bU) 

where 

n  =  Y(q+Y)_1x  +  qt((q+Y)(q+Y+kt))_1  x  U  +  q(q+Y+kt)  m  1 
a  =  Yq(q  +  r)"1 
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b  =  q2t(q+Y)_1 (q+y+kt)'1 

X  ~  N(ml,(q  +  y) I  +  tU) 

Lemma  2 . 5 .1 .  Let  Y  ~  N(|i  +  pl_,aT  +  bU)  wi  th  m  f  IR  k ,  p  C  IR  ,  a  0  and 

b  •  -  a/k.  Then  there  exists  a  random  vector  Z  -  N(p,al)  such  that 

k  k 

h(Y)  -  h(Z)  everywhere  for  every  translation-invariant  n  :  IR  >  ]R  . 
Proof.  (See  Miescke  (1979)). 


With  this  lemma,  it  is  easy  to  get 

P^x)  -  P(ui  =  n[-k]|x) 


(u) 


where  is  the  normal  distribution  with  mean  y  and  variance-covari- 

ance  matrix  V. 

We  can  rewrite  p^x)  as 


Pi  (x)  =  /  ^  +  (qTq+7T)L(xi  '  xj))d*(t)  • 

2  2 

I et  1  Hq,  q  =  o  /n ,  we  have 


P:(x)  =  /  ll  <t(t  +  ( 

j/i  ,  ol  +  2, 

n  (  n  °0' 


0 

2  2 


) ** ( x i  -  Xj))df(t). 


The  above  expression  is  exactly  the  same  as  that  of  the  independent  prior 
Case  I,  lb. 


Case  V.  Under  normal  assumption  as  before,  but  suppose  -.’s  are  unknown 
and  Thai  neither  .^'s  nor  n-'s  are  all  equal. 


Suppose  we  have  no  prior  information  about  (;.  ,o ) ,  for  each  mdividu- 
al  population  ^  assign  prior  p(u.,o-)ro^  then  we  have  (See  Bo*  and 
Tiao  (197  3))  that  the  posterior  density  of  im  ,  given  X  r  (x .  j , .  . .  ,  x 


P(ui |xi ) 


(s  -//FT.")” 

_ 1  _ 

B  ( i,  ,  i )  J'o . 


n .  (u  •  -x  • )  *1 ) 

11  1 


where  s^  is  the  sample  variance,  B(*,*)  is  a  Beta  function  and 


Hence 


Mi'xi 

p(t,  =  —  --1-  j  x  . )  --  — - 
s/  *ni  B(.' 


B  ( .  v  ■  ,  )  /v  1- 


t,  I ) 

(1  +  '  )  ’ 


which  is  the  density  of  tne  student's  t  distribution  with  ,o  (  n.  -  1) 
degrees  of  freedom. 

Using  this  result  we  can  write  the  formula  of  p^(x)  by 
P i  ( x )  -  P((<i  ■  Up  V  j/1  ;x) 

s./t'n.  x.-x. 

/  :  1  Tv  •  ( t  '  . 1  +  J  )dTv  (t) 

jO  Sj/Vnj  y/..j 


where 


; .  =  n .  -  1 ,  i  =  1 , . . .  ,k 

li 


’  .  s  .  =  j  (  X  •  -  X  •  ) 

ii  L  i  i  r  i 
r  -1 


T„.  is  the  c.d.f.  of  t  distribution  with  .  degrees  of  freedom, 
l  i 

When  v/s  are  large,  t  distribution  approaches  normal  distribution,  home, 
fnr  large  n.,  i  1 . k,  we  can  replace  T  by  !. 
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Case  VI.  Suppose  we  are  interested  in  finding  a  subset  which  contains  the 
population  with  the  smallest  variance;  i.e.,  we  define  the  best  population 
as  the  one  with  the  smallest  variance,  and  suppose  that  we  have  no  prior 
information  about  o.  In  this  case,  it  is  reasonable  to  assume  that 

p(;i ,a  )..o  ^ ,  i f  ,1  is  unknown 

p(.i ).. .1  \  if  ;i  is  known . 


I.et 


9  ,  9 

n;,  v-S^ :  =  }  (X;r-.,)  if  ..  is  known 

k 


J'  'J"J  (,L; -|  '  Jr 


■i.  -  n  .  -  1 ,  v.S2  -  i  (X  .  -  X  .  )2  if  p.  is  unknown ,  n  .  •  1  j  =  1 ,k 
.1  .1  j  J  r-]  Jr  J  J 

S?  ■  . Sk>.  X  -  <X,1 . Xln, . Xkn  > 

i  k 

2  2 

and  Y  be  the  random  variable  with  c.d.f.  ;<  which  is  the  v  distribution 
with  degrees  of  freedom. 

Then  for  either  case  (u  known  or  unknown),  we  have 


P i  ( x )  =  p('^  "  l'[  |  ] !  x  =  x) 


=  P(.'j  •  Oy  X  j/i  IX  =  X) 


,S2  ,.s2  v.S2 

=-  p(  y  •  v  (-J-J 

V-S". 

J  1  11 


-'j  (  -J-  -2),  V  jj*i  ;S2  -  s') 


P(Y 


^2 

Y  (  J),  \  j/i  |S2  -  s?) 


2 

v,  s  . 


f  (u  'j4)d-*  (u) 

U  .|Y  i  v  jSi  '  i 


J 


Ill 


-  /  •;  /(u  )  f  n  ,  ^ .  n  » i . 

o  j/i  •’  -  1  k 

With  these  p-j  (x) , . . .  ,P|.  (x)  we  can  apply  Bayes-P  i~ules  .  iwut  t,r,:!v. 

B  B  * 

Lemma  2.5.2.  In  Case  VI,  ^  and  ,,f,R  are  just  a.e.  and  (stale]  tr.insl.. 
tion  invariant. 


*  Here  the  definition  of  the  "just"  property  for  a  selection  rule  ■ 


Ms2’)  if  S? 


.2 ' 


s^ 


\i 


B  M 

2.6.  Comparison  of  Selection  Rules  and  in  the  Normal  !  o<  at  .un 
Parameter  Ca se 

We  have  k  normal  populations  with  a  common  known  variam >■  < 

common  sample  size  n.  For  this  case  Gupta  (1966)  proposed  uni  '  .n 
procedure 

Select  a.  iff  X.  •>  Xr.  ,  -  d  i  •  l,...,k  when  ' 

1  1  /n 

is  to  be  determined  by 


inf  P(CSjvM)  P* 

"  ( 

and  l  is  the  parameter  space. 

M  ★ 

We  will  show  that  o  (  (  , P  )  where  r  is  the  locally  unito- 

* 

prior  distribution,  for  fixed  D  and  k,  let  <1  be  determined  by 


/  i-k_1(t  +  d)dt(t)  -  P*  . 


Let 


fall  possible  observed  values, 


IX  C,  !x[k]  -  d  ;■  •  Xl1] 

v  n 


■i  ix  c  ;c  1  x [ i - 1 -j  •  x : k ]  ■ d  7  x[i j f ’  ?  1  - k 


(X  (  !X,  ,  ,  =  Xr  •  ■  !  ■  Xr  -  d  -  X,  ;  c 

I  L  U  L  1  _  1  J  l  k  J  /_  L  1  J 

v  n 


iXM]  ‘  xh-i]  ‘  x[^j  •  d  /n  =  x[ij  ^  x[i-i],c-  - 


then  we  have  the  following  theorem. 


■k  1  * 


theorem  li.6.1.  Given  a  number  P  (  ■  P  ■  I)  and  locally  uniform  prior 


for  each  population  X  ~  x  (  X  . ,  then 


P(CS..'\  X  -  x)  ■  q  (i) 


*  L  _  i  *  * 

q  (D=  (1  -  P  )  +  P  ■ 


/V,  N8r,.p*) 


I’rnol.  It  is  sufficient  to  show  that 


inf  p  •  (x )  -  q  ( ' )  -  k_-|  (i  -  P  )  +  P  . 

x  (  .*c  .  •  -  i  ^ ' J 

i 


,\\i  int  P(rs 


1  Of', 


=  1  -  (i  -  1)  /  :■{ t-d )  Jk"2(t)d*(t) 

—  rxt 


=  (k  -  i)  /"  ok_2(t)<f(t  -  d)d-:(t) 


+  /  ;k_1  (t  +  d)ck’(t)  (2.6.3) 


(2) 

The  superimum  of  (2.6.3)  occurs  when  x  6^,  :  .  The  last  equality 

follows  from  the  identity 

(k  -  1)  /  ck-2( t)<:>(t  -  d)df-(t) 

=  1  -  /  tk-1 (t  +  d)dt(t) , 

which  can  be  shown  by  the  integration  by  parts.  By  (2.6.1),  the  second 
term  of  (2.6.3)  equals  P*;  then  use  the  integration  by  parts  to  the 
first  term  of  (2.6.2),  we  get 

inf  l  Pm(x)  -  k'4  [1  -  P*]  +  P*  (2.6.4) 

x  J.=i  L  J 

=  q*(i )  . 


M 

Remark  2 .6.1.  If  the  procedure  «;  selects  ^  ^  only,  i.e.  X  =  xty 

B  B 

then  by  Theorem  2.6.1  we  have  p^-j(x)  ^  P*  so  that  or  selects 


(k ) 


only.  But  the  converse  is  not  necessarily  true. 


B  M 

Remark  2.6.2.  For  the  case  k  =  2,  =  il»  a.e.  For  any  given  X  -  x; 

M  B 

il  x  (  then  Pj-pj(x)  s  P*,  hence  v  and  tjj^  select  the  population 

m 

associated  x^j.  If  x  C  ,  and  x^j  -  d  -  :  Xp-j  then  •  and 


B 

VNR 


select  both  populations 


P(X 


[2]  ' 


and  >9 


Since 


X[l] 


)  =  0, 


we  have  y 


B 

NR 


M 

v  a .  e .  . 


Remark  2.6.3.  The  above  Theorem  and  Remark  2.6.1  gives  us  a  lowe 

k 

bound  on  the  value  of  ]_  Prr-i(x),  over  all  x  t .  The  exact  v.il 

c=  i  -*  1 

of  the  difference  of  the  selected  sizes  between  ^  and  depend', 
the  observations. 


u  .-a 

2.7.  Applications  to  Select  max  •>,  ,  o.  =  .  f0r  Normal 

1  <i  tk  1  ’j 

O 

Distribution  N(u.,  at),  i  -  l,...,k 

Let  be  k  independent  normal  populations  with  mean  e 

2 

variance  o . ,  both  u.  and  o.  are  unknown.  For  the  goal  of  finding 
random  subset  wtiich  contains  the  population  with  maximum  a. 

i 

for  some  given  constant  a,  we  assume  that  apriori  (i..,  .  ),  i  I 

are  independent.  Suppose  we  have  n •  independent  observations 
X^,...,X.fi^  from  ,  and  let  X^  bo  their  sample  mean,  i  1,..., 

Let  Yj,...,Y  be  i.i.d.  -  N(„,  ).  If  no  prior  infonn.it  ion 

available  to  r),  we  could  assion  a  locally  uniform  prior 
p  ( ;  i ,  ')  .<  •  ^  to  o),  (see  Box  and  Tiao  ( 1 97  ft ) ) .  And  the  post. 


r 


IIC 


r  »r  i 


.  «md 

i 


,) 


_ k 


I-  . 


•r  i<>» 


joint  distribution  of  ,i 1  :  n  -  a  and  n,  qivon  obsrrval  ions 
Y  r  y  ■  (y | . yn)  is  given  by 


nn 


P(w\  ~'\y)  =  kn  ^n+1'1  exp  { - --x-  |\>s2  +  n(y‘  -  u ' ) 2 ]  f 

2  c; 

where 


y‘  =  y  -  a,  y  =  7  y./n 

=  1  1 

vs2  =  l  (yi  -  y)2,  V  n  -  1  ^  (2.7.1  ) 


k  =  /2"  ['2  "  ( p )  ] " 1  (  2"  ) 2  ■ 

Let.  t  -  /n  (u-a)/n,  with  (2.7.1)  the  posterior  distribution  of  ?,  given 
Y  -  y  is 

i»  ( '  l Y  =  y)  =  p(f.|t) 

v  -i  v 

=  f2?  r(pr](  - -)2  exp  {-  (  ^%-)t  '(v)I  ,(  ^-v/o) 

2  v+t2  2  v+t2  v_1  (,+t‘r/? 


where 


t  /n(y  -  a)/s,  v  n  -  1 


I  (x)  -  /  ( >2”  r(v))  1 u exp  {-  l  (u  +  x )2  idu . 

0  2 


Now,  let  p(;w,  '•)  ■  be  the  assiqnerl  locally  uniform  prior 

to  (i-j,  .).  Then  let  x  =  ( x ,  ( - ,x]n  . x^  ),  we  have 

1  'k 

i . 


Pi  (x)  --  P ( -■* i 


[k] 


) 


i  •  “  i  <  •  -  d 

-  P(  1  max  (  J  )  lx) 

’  i  1  •  j  ■  k  '  i 


=  P( 


. 1 


n . 


i 


f, . 
J 


V  j  /  i  1 1 ) 


=  /  •!  G.  (/ --i  z  |  t)d  Gf  (z  1 1) 
j/i  '  j  i  '  i 

=  /  II  G_  (z!t)d  G.  (z|t)  if  n.  =  ...  -  n,  n, 
j/i  j  i 

where  Gf  is  the  posterior  c.d.f.  of  r,  given  x  or  t. 

By  (2.7.2),  the  Bayes-P*  procedure  is  completely  specified. 

If  the  prior  distribution  for  (;>,  <■.)  is  the  conjugate  distribution 
(see  Raiffa  and  Schlaifer  (I960)),  then 

p(u,  0  ,x  exp  {-  n 1  (s-m1  )*'  •  ^  •  exp  I-  1  ?  {  ’  / 1 ) 

2  <  t  ‘  7 ' 

-  p(u|n)  p(e) 

that  is 

p ( .■  | a)  -  N ( m ‘  ,  ,/./n  '  )  ,  n'  l.) 

nl  S  vV  2 

p(o)  -  -~~2  ■  -v.  »  V  ,  v  0. 

Let 

,  nx  +  n'm'  .,  . 

x  =  —j —  ,  x  is  the  sample  mean. 

u2  -  1 (n  -  l)s?  +  v ' v 1  +  [(nn')/(n  +  n')J(x  -  m')fi/  * 

V*  r  (n  -  1 )  +  V 1  +  1 

•*  -  (n  •»  n')1/2(,  -  a)/- 
t*  -  (n  +  n'  )1/,?(x‘  -  a ) / u , 

the  posterior  distribution  of  r*,  given  x  is  p  ( ■’  * '  x )  ■  p(  *'t‘)  which 
has  the  same  form  as  p(r.!t),  but  replace  •  ,  t,  ■  by  * ,  t*,  *. 


I  hu'.  tor  the  i  oil  I  n<|'i  t  f  I'l'ini'  t  wi1  net 


'/  1 


i  ^  ,  i  t  *  1  j  *  n  1 

,t/i  i 

* 

where  (i  t  is  f  he  posterior  f  .  !.*.  of  <jiven  x  or  t. 

Note  that  (i’,/.  t)  has  the  same  form  as  (2  ./.?)■,  but  I »  <• 

,  t  by  *,  t*. 

A(>|>  1  it  at  inns  to  Poisson  Pis  tributions  an  <  t’o_i->0"  <•»•,,,,  .  ■■ 

I.  Poisson  Distributions  Case 

'•oppose  that  ^  a  no  k  independen  t  Poisson  oopu  Kit  i  <>ns  ,  w1'-" 

t  ho  independent  observations  X. j .  ,X.  from  -  have  the  po issun 

i 

•tensity  with  parameter  •  ;  denoted  by  P(*  ■•).  i 

lot.  V| . Y()  be  i.i.d.  with  p(*  ■)■  It  we  use  non-intoniaf  i  .<■ 

-  )  ft 

prior  p(  •  )  .  '■  '  ‘  (Box  and  Tiao  (i°7‘)).  then  >ii  ven  v  y  ;  t 

we  hrive  the  posterior  density  as  follows: 


P 


•y 


nv- 

r  ■ 


exp  ( - n ■ 1 


where 


y 


n 

i  1 


y.  and 


•  r  v 


■> 


We  see  that  2n  \y  '-2ny+l  ’  the  ch’~sc)udro  distribution  with 
2ny+l  degrees  of  freedom.  Hence  by  using  rion-informati ve  prior 

p(\  . )  '■  for  each  population  -  we  have 

pi(x)  p(S  =  '[k]  x) 

=  1  V.  (2  -  )d  (7) 

-•  j/i  J  j  1 

where 

n  ■ 

1 

■  "  2 n  .  x  +  1  ,  x  .  =  }  x  .  •/ n  . . 

1  it  1  •  l_ ,  1  j  1 

o  ' 

If  n,  .  .  .  -  n.  ,  then 
I  V 

p5(x)  =  /  ■■  ■/:  (z)d  x?  (z). 

0  j/i  7  i 

With  p.(x),  i  we  can  apply  Bayes-P*  selection  rule 

.  ^  and  easily  to  select  a  subse*:  which  contains  the  popul.it  : 
with  the  largest  parameter  *  .  On  the  other  hand,  if  we  arc  inter 
in  selecting  the  population  with  the  smallest  parameter  •  ,  t hi  n 

p.  (x)  /  •  [1  -  /  (z  "jHd  /  { ? ) 

0  j  i  i  "  i  i 

■  /  "  [  1  -  ( z )  _ld  ■ '  ( Z  )  i  f  n  ,  .  .  .  n 

0  .1  /  i  j  i 

In  this  case,  the  simulation  results  tor  selection  proced.ir. 

arid  J*  are  tabulated  on  Table  VII. 

W< 
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2.8.2.  Poisson  Processes  Case 

Suppose  we  have  k  independent  Poisson  processes 

IX^(t)l,...,|X^(t)!  with  expected  arrival  times  equal  to 

'  .  -  ,  respectively.  Hence  for  the  processes  ;X^(t)l,  the 

'l  k 

probability  that  there  are  in.  arrivals  until  time  t.  is 


p(X(l)(ti)  -  mi|Ai,  t.) 


m. 

(t.\.)  1  -t->. 


l  l 

in. ! 

l 


i  i 


i  =  1,  2 - ,k,  ni .  =  0,  1  ,  2  , .  .  . 

If  there  exists  no  prior  information ,  then  we  use  the  non-i nformati ve 
-1/2 

prior  p(\.)  -  a ■  for  all  processes.  Therefore,  we  get  the  posterior 
density  function  of  .U  ,  given  (m. ,  t^)  as  follows: 

p(Ai !X(l ^(ti )  m. ,  ti )  p ( x  .  im. ,  t . ) 

m.  +  1  -1 

( t .  •  . )  1  ‘  - 1  •  A  • 

=  ’--1  t.  e  1  \ 

“'(’"i  +  ?} 


Ihus  8fj‘-  has  •'  distribution  with  2nu  H  degrees  o *’  freedom,  given 

l  iie  number  m.  of  arrivals  before  tine  t. . . 

l  i 

let  hi  •  (in  |  , .  .  .  ,111^ )  and  t  (t^ . 1  ^ ) ,  then  it  <  an  tie  shown  th.it 

the  Poisson  process  <'X^^(t'  has  the  maximum  parameter  (or  minimum 


mean  waiting  time),  given  (m,  t.)  is 


Here  we  list  two  special  cases  which  are  of  interest. 

(a)  Observations  of  all  processes  are  obtained  in  a  common  t  He 
interval  [s.,  t  +  s .  Since  Poisson  process  is  stationary,  we 
can  assume  that  s^  =  0,  and  t^  =...=  t^  =  t.  In  this  case 

Pi('"’  t}  =  /  *L.  +  l(y)d  4m.+1(y' 

0  1  i 

which  is  independent  of  t. 

(b)  All  fir's  are  equal,  i.e.  we  fix  m  first,  then  net.  observation-. 
Hence 

Pi  (,n*  ^  =  v2m+l^y  t^'1  ,;?n)+l  •  ’ 

There  is  an  alternative  way  to  approach  the  cases  (a)  and  (hi. 
Let  T .  be  the  waitinq  time  of  the  nth  arrival  in  the  ith  process, 
then  T..  has  a  gamma  distribution  with  density  given  by 

n .  - 1  -  \ .  t 

Pit)  •=  .(.;,;v  (y)  1  e  1  t  -  0. 

I/O 

If  we  have  only  nbn-informati ve  prior  p(\)  '  ;  then,  given  in.  ,u 

? 

t  -  ,  ('th.  nas  posterior  distnputian  ..,p)  +1 ,  therefore  the  i.-r  -ila 

of  p^(m,  t]  we  net  here  is  exartly  *  he 


same  as  before1. 


Remark  P.B.l.  Under  non- inf onnati  ve  prior,  in  comparing  the  subse* 
selection  problem  in  k  Poisson  distributions  with  the  problem  in  k 
Poisson  processes,  it  is  easily  seen  that  Poisson  distributions  mode  I 
is  a  special  case  of  Poisson  processes  model,  namely,  t..  =  n.  where 
n.  denotes  the  sample  size  of  the  ith  Poisson  population. 

P.R.3.  Relation  Between  Selection  from  Poisson  Processes  and 
Selection  from  Populations  with  Gamma  or  Exponential 
Distribution 

Suppose  we  have  k  independent  populations,  the  ith  population 
havinq  the  qamnia  distribution  with  parameters  .  =  m.  (known), 

.  “  1/-.  (unknown).  Since  the  random  variable  1\  ,  the  waiting  tine 
until  m.  arrivals  in  a  Poisson  process  with,  parameter  has  a 
qamma  distribution  with  parameters  .  =  m.  ,  -!/<..  If  tne  m .  1  s  are 

qiven  and  if  the  goals  for  both  selection  problems  are  the  same, 
namely,  to  select  a  subset  containing  the  population  (process)  with 
the  largest  parameter  »,  then  it  is  easily  seen  that  these  are 
identical  problems.  Note  that  in  the  selection  problem  of  Poisson 
processes,  nr’s  might  not  fy>  the  proassiqned  .  iuc-s  but  are  given 
random  observations  whenever  t . 's  are  proassiqned  values.  In  this 
1  ii'.ig  the  selection  problem  of  Poisson  processes  is  different  fro.,, 

I  he  *  of  the  gamma  distributions. 

If  the  process  asso,.  iaf  e-:  with  the-  minimum  parameter  ■  (or  the 
mi-imum  waiting  time)  i t  he  »  ,  Mien  the  posterior  probability 

l  ;  \ 

of  iil'e,  esc  X '  { 1 1-  t  a  i-*  t  ’ 


!n  ‘  m  antilogous  to  f  he  one  off  a  i  ne< 1 


before  with  the  modifications  that  the  integrand  function 


j/i 


2 

‘2ntj+l 


<y  tJ) 

1 


of  (2.8.1)  is  replaced  by 


t'  -  4,0  tJ  n. 

j/i  i 

2.9.  Comparison  of  the  Performance  of  Cw,  and  ;Ml‘J 

Let  n  .  i  -  l,...,k  be  k  independent  populations,  where  .  h.i  .  i.hc 
associated  c.d.f.  F(x,  9.)  -  F(>;  -  0.)  with  unknown  loca  t  i  on  |i<)  i  .  him  *  t  *  *: 
‘K  .  Let  f (x,  ■'.)  =  f ( x  -  :'.)  be  the  p.d.f.  The  goal  is  to  ’in.' 
small  (nontrival)  subset  which  contains  the  best. 

The  following  subset  selection  procedure  based  on  vi.vli 
medians  is  due  to  Gupta  and  Singh  (1980). 

MFD 

il>  :  Select  ■■■  ■  if  and  only  if  X..  Xj-^j  -  d 

where  X.  is  the  median  of  the  2m+l  random  observations  from  popu  1 t  i . n 

it.  and  X  =  max  X..  The  value  d  is  determined  hv  the  following 

[k  ]  ■} 

equation  so  that  the  P+-condition  is  met. 


where 


/  G(u  +  ri) '  u(u)du 


P* 


r)  ( u ) 


(?m+l ) ! 
(m!)? 


[ r ( o ) ]  ’r '  -  r(u) ]nif(u) 


I 


G(u)  =  Ip(uj(i!i  +  1  ,  1:1  +  1) 

I  (p,  q)  is  the  incomplete  beta  function. 

In  this  section  we  use  Monte  Carlo  simulation  techniques  to  compare 

B  B  M  MED  . 

the  performance  of  selection  procedures  .■  ,  •.  and  .  in  the 

M  f-i  y  Q 

normal  means  problem.  Because  both  rules  -  and  .  are  not  based 

on  any  prior  information  about  trie  unknown  parameters,  we  assume  that 

the  prior  distribution  i  for  botn  and  is  locally  uniformly 

M 

distributed.  Since  the  selection  procedure  '  satisfies  both  the 

P*-condition  and  the  poster ior-P*  condition  wrt  the  locally  uniform. 

B 

priors,  it  makes  sense  to  compare  the  Bayes-P*  procedures  .  and 
i.;  m  m  Mr  n 

mj,  with  .  and  compare  .  with  .  '  in  conns  of  efficiency  whicn 

is  the  ratio  of  the  probability  of  a  correct  selection  to  the  expected 
selected  size,  for  studying  the  robustness  of  these  four  rules, 

p,  1.1  MCI) 

'  and  .  ,  we  change  the  true  distribution  to  non-normal  dis- 

trituitions,  namely,  the  logistic,  Laplace  (the  double  exponential) 

•>nd  the  gross  error  model  (the  contaminated  distribution),  but  keep 
the  selection  procedure  unchanged  (i.e.  still  baseo  on  the  normal 
as  .-imp!  ion) .  Inc  Mon  t  e  Carlo  simulation  results  for  both  equal  dis- 
tm.es  of  the  parameters  and  slippage  cases  are  tahul  a  tod.  In  the 
.  i  rn  1 1  a  t  i  on  stud?  all  generate-:  r.injnm  variables  arc  adjuste.*  to  ha  v 
-.armin' u  1.  i  a<  h  t  t  :•.(*  w*  tjpr.f  cat  °  ‘i.e  random  variables  w’ t  h  tit1 
'ir/i'M  d  i  s  t  r  i  bu  t .i  on  of  cacti  p»s;a  In  '  .('■  ,  ‘hen  apply  f  h»  seie-  »  !<••••  pro- 


■  •  d-e'es .  The  simulation  process  i i-  repeated  1  DO  •  i  o  4 ■  *»• 
/.is  I  able  .  The  relative  !  reu  icnc  .  o'  se  I  is  t  i  ir;  the  pnp.ilat 


<M  -  *1  r  "i  ts 


used  as  an  approximation  to  the  probability  of  selectin'!  the  om'U 
lation  i  * .  The  sum  of  relative  frequency  of  selectin')  each  pup.i 
lation  "j,  i  -  is  treated  as  an  approximation  of  the  '■«■)» 

selected  size.  The  efficiency  EFT  of  each  selection  pmcnlure  i  . 
approximated  by  the  ratio  of  relative  frequency  of  selectin'!  the 
best  one  to  the  expected  size.  The  simulation  results  indim'e 
that  in  all  cases  we  have  the  performance 

5  B  M 

V  r  ■  NR  >  1>  • 

It  should  be  noted  that  in  the  above  comparison  cf  the  (.'t-r ! om.im  i 

we  restrict  attention  to  these  rules  which  satisfy  the  pos' e-'int  -! 

condition,  for  small  sample  size,  the  efficiency  cf  rule  ,  ten: 
MEO 

to  be  larqer  than  ,  under  P*-conriition. 

Remark  2.9.1.  The  Laplace  distribution  has  the  density  f unct 

( x  -  ■  )  -  -  X  • 

for  which  the  variance  is  2. 

The  loqistic  distribution  has  the  density  fniu  t ion 

'•(>•  -  •  ■  =  " 

'>e  '  ')' 

for  which  the  variance  Var  ;v)  -  t 

The  press  error  model  we  ns»ei  hat  'hr-  density  fan-  ■  ion 
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i  he  efficiency  o'  a  selection  procedure  .  i  . !«■  i  i r. f*>.l  t •  y 


r  r 


I1  (rs‘ . } 

:  (s'.) 


where  [.  ( is  the  expectm  selected  size. 


Discussion  and  Conclusion 

for  Table  VII  1.1  and  TabiC  VI II...'  (equal  distances  case;  the 
is  .99  and  .90  respectively,  tne  common  sample  size  n  -  0 .  k  s. 
It  the  k  populations  have  normal  Msf  ritv.it  ions  with  the  unknown 
parameter  configuration  +(k-l;'}.  com  on  variance  1.  frc 

both  tables  the  performance  eased  on  either  the  efficiency  or  the 
experteit  selected  size  is 


i1  t  tie  nosterior-P*  condition  is  considered,  and 
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under  t  he  P*-i  nmti  t  i  on  . 
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have  normal  distributions  with  unknown  parameter  conf i pur  a n .a--, 

(  j  ,  common  variance  1.  From  both  tables  the  pert orman- 

based  on  either  the  efficiency  or  the  expected  selected  size  i  . 

B  B  M 

v  >  >;r  >  v 

if  the  posterior-P*  condition  is  considered,  and  if  ->n  1 

>  , 

under  the  P*-condition. 

Note  that  in  both  equal  distances  and  slippaoe  cases  when 
that  is  the  population  means  are  not  very  close,  the  procedure. 

g 

i'NR’  wrt  the  locally  uniform  priors,  always  satisfy  not  only  tin 

posterior-P*  condition  but  also  P  (CSt.B  or  ,B  )  P*,  and  the  ■ 

11 

ed  selected  size  of  the  selection  procedure  ,B  or  in  mm  h  1. 

NR 

than  the  selection  procedures  ;M  and  .,MED.  For  example,  in  *  hi- 

★ 

ma 1  equal  distances  case,  P  .99,  k  =  5,  ten  -  4. 

C(<UMED)  -  E(SUjjR)  ^  P.38P; 

★ 

in  the  normal  slippage  case,  P  =  .99,  k  =  5,  * /n  4 
E(S|  :MED)  -  t(S'-.BR)  -  1.560. 
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